HOMOGENIZATION OF THE NEUMANN PROBLEM FOR 
ELLIPTIC SYSTEMS WITH PERIODIC COEFFICIENTS 



T. A. SUSLINA 

Abstract. Let O C R'' be a bounded domain with the boundary of 
class C^'^. In L2{0;C"), a matrix elhptic second order differential 
operator AN,e with the Neumann boundary condition is considered. 
Here £ > is a small parameter, the coefficients of An.e are periodic 
and depend on x/e. There are no regularity assumptions on the co- 
efficients. It is shown that the resolvent {An,£ + converges in 
the 1/2(0; C")-operator norm to the resolvent of the effective operator 
A'if with constant coefficients, as e — > 0. A sharp order error estimate 
\\{An,e + ^I)^^ — iA% + \I)^^\\l2^L2 ^ Ce is obtained. Approximation 
for the resolvent (.4jv,e + A/)~^ in the norm of operators acting from 
LiiO; C") to the Sobolev space H\0; C") with an error 0{y/e) is found. 
Approximation is given by the sum of the operator {A% + and 
the first order corrector. In a strictly interior subdomain O' a similar 
approximation with an error 0(e) is obtained. 



Introduction 

The paper concerns homogenization theory for periodic differential oper- 
ators (DCs). A broad literature is devoted to homogenization problems in 
the small period limit. First of all, we mention the books [BeLPa], [BaPan], 
[ZhKO]. 

0.1. Operator-theoretic approach to homogenization problems. In 

a series of papers [BSul-4] by M. Sh. Birman and T. A. Suslina an operator- 
theoretic (spectral) approach to homogenization problems was suggested 
and developed. By this approach, the so-called operator error estimates in 
homogenization problems for elliptic DCs were obtained. Matrix strongly 
elliptic DCs acting in L2{^'^; C") and admitting a factorization of the form 

Ae = b(Dyg{x/e)b(D), e > 0, (0.1) 

were considered. Here g{x) is a periodic bounded and positive definite 
matrix-valued function, and b(D) is a first order DO. The precise assump- 
tions on ^(x) and b(D) are given below in §1. 

In [BSul-4], the equation AeUe + u^ = F with F e L2{W^; C") was consid- 
ered. The behavior of the solution for small e was studied. The solution 
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Us converges in L2(M'^; C") to the solution uq of the "homogenized" equation 
A^uo + uo = F, as £ ->■ 0. Here = b(D)*g^b(D) is the effective operator 
with the constant effective matrix ^f'^. In [BSul,2], it was proved that 

In operator terms, it means that the resolvent {Ae + 1)^^ converges in the 
L2{W^; C")-operator norm to the resolvent of the effective operator, as e ^ 0, 
and we have 

II (A + I)-' - + I)-'\\L,iR'i)^L,(R'') ^ Ce. (0.2) 

In [BSu3], more accurate approximation (including the corrector) for the 
resolvent {Ae + in the L2(M°'; C'^)-operator norm with an error term 
0{e^) was obtained. (Here we do not go in details.) 

In [BSu4], approximation of the resolvent of Ae in the norm of operators 
acting from L2(M'^;C") to the Sobolev space H^{W^;€V-) was found: 

||(^, + /)-! - (^0 + /)-i - eK{e)\\L,(^<i^^H^i^^'>^) ^ Ce- (0.3) 

this corresponds to approximation of in the energy norm. Here K{£) is 
a corrector. It contains rapidly oscillating factors, and so depends on e. 

Estimates of the form (0.2), (0.3) are called operator error estimates. 
They are order-sharp; the constants in estimates are controlled explicitly in 
terms of the problem data. The method of [BSul-4] is based on the scal- 
ing transformation, the Floquet-Bloch theory and the analytic perturbation 
theory. 

0.2. A different approach to operator error estimates in homogeniza- 

tion problems was suggested by V. V. Zhikov. In [Zhl, Zh2, ZhPas, Pas], 
the scalar elliptic operator — div5f(x/£)V (where ^'(x) is a matrix with real 
entries) and the system of elasticity theory were studied. Estimates of the 
form (0.2) and (0.3) for the corresponding problems in were obtained. 
The method was based on analysis of the first order approximation to the 
solution and introduction of an additional parameter. Besides the prob- 
lems in M*^, homogenization problems in a bounded domain O dW^ with 
the Dirichlet or Neumann boundary condition were studied. Approximation 
of the solutions in H^{0) was deduced from the corresponding results in 
W^. Due to the boundary influence, estimates in a bounded domain become 
worse and the error term is 0(£^/^): 

\\A^^,-{Al)-^-eK,{e)\\L,i^o)^HHo)^Ce'/\ b = D,N. (0.4) 

Here AD,e and AN,e are operators with the Dirichlet or Neumann bound- 
ary conditions, A%, A% and Kd{£), KN{e) are the corresponding effective 
operators and correctors, respectively. 

In [ZhPas], an estimate \\A^^ — {A^)~^\\l2^L2 ^ Ce^^"^ was obtained as a 
(rough) consequence of (0.4). Refinement of this estimate is a natural prob- 
lem. In [ZhPas], for the case of the scalar elliptic operator — div5(x/£)V 
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(where g{x) is a matrix with real entries) with the Dirichlet boundary con- 
dition an estimate for H^^^^ — (^^j)"^ ||l2->L2 order £'2d-2 for d ^ 3 and 
of order e\ loge| for d = 2 was obtained. The proof essentially relies on the 
maximum principle which is specific for scalar elliptic equations. 

Similar results for the operator —divg^x/e)^/ in a bounded domain with 
the Dirichlet or Neumann boundary conditions were obtained in the papers 
[Grl, Gr2] by G. Griso by the "unfolding" method. In [Grl] estimate of the 
form (0.4) was obtained, and in [Gr2] the sharp order estimate 

WA'e - K)"' KiO)^L,iO) ^Ce, b = D, N, (0.5) 

was proved for the first time (for the same scalar elliptic operator). 

0.3. Operator error estimates for matrix elliptic operators in a 
bounded domain. In the recent papers [PSul,2], [Sul,2], the Dirichlet 
problem for the equation Ad^sU^ = F with F G L2{0;C"-) in a bounded 
domain O C with the boundary of class C^'^ was studied. Here AD,e 
is a matrix operator of the form (0.1) with the Dirichlet condition on dO. 
In [PSul,2], estimate of the form (0.4) for the operator AD,e was obtained. 
The method was based on using estimates (0.2), (0.3) for the problem in 
and on estimates for the discrepancy (of the boundary layer type) ; some 
technical tricks were borrowed from [ZhPas]. 

In [Sul,2], the author succeeded in proving a sharp order error estimate 
of the form (0.5) for the same matrix operator Ad,s- The problem was to 
obtain an estimate of order e for the L2-norm of the discrepancy; for this 
purpose, estimate of the form (0.4) was applied and the duality arguments 
were used. 

In a recent paper [KcLiS] , homogenization problems for uniformly elliptic 
systems in a bounded domain with the Dirichlet or Neumann boundary 
conditions have been studied. Under the assumptions that the coefficients 
are real-valued and Holder continuous the authors obtained estimate of the 
form (0.5) in the corresponding problems. 

Note that the class of strongly elliptic operators (0.1) that we consider is 
wider than the class of operators studied in [KeLiS]. Moreover, we do not 
impose any regularity assumptions on the coefficients, which extends the 
of possible applications. 

0.4. Main results. In the present paper, analogs of the results from 
[PSul,2], [Sul,2] for the Neumann problem are obtained. Note that the 
homogenization problem with the Neumann boundary condition is more 
complicated than that with the Dirichlet condition. Main difficulty is re- 
lated to the fact that the boundary conditions in the initial and homogenized 
problems are different. Indeed, the conormal derivative of the solution con- 
tains rapidly oscillating coefficients, while in the homogenized problem the 
conormal derivative has constant coefficients. (For the Dirichlet problem, 
boundary conditions in the initial and homogenized problems are one and 
the same.) 
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We study matrix DO's An,s in a bounded domain O G with the 
boundary of class C^'^. The operator An,£ is given by the differential ex- 
pression (0.1) with the Neumann condition on dO. The effective operator 
is given by the expression b(D)*g^b(D) with the Neumann condition. 
The behavior of the solutions of the equation {An,£ + A/)u£ = F for small 
£ is studied. Here F G L2{0; C"-), and A is subject to the restriction which 
ensures that the operator AN,e + A/ is positive definite. The case where 
A = 0, which is important for applications, is studied separately in §9. In 
this case, additional orthogonality conditions on F and are imposed. 

In operator terms, the following estimates are obtained: 

||(^^,, + XI)-' - {A% + Xir^\\L,iO)-.L,iO) ^ Ce, (0.6) 

||(^^,, + XI)-' - {A% + XI)-' - sKM{s)\\L,io)^HHo) ^ C-eV^. (0.7) 

Here Kj^is) is the corresponding corrector. The form of the corrector de- 
pends on the properties of the periodic solution A(x) of the auxiliary prob- 
lem (1.7). In general case, the corrector contains an auxiliary smoothing 
operator. If A is bounded, the corrector has a standard form. Besides ap- 
proximation of the solution in H'{0; C"), we also obtain approximation 
of the flux Pe = g'^b{'D)u^ in L2(C';C™'). For a strictly interior subdomain 
O', the following sharp order estimate is proved: 

\\iAN,e + Xiy' - {A% + XI)-' - eKN{e)\\L,iO)^miO') ^ Ce. (0.8) 
The author considers estimate (0.6) as the main achievement of the paper. 

0.5. The method is based on using estimates (0.2), (0.3) for homogeniza- 
tion problem in and taking the boundary influence into account. Main 
difficulties are related to estimation of the " boundary layer correction term" 
We which is solution of the Neumann problem for the homogeneous equa- 
tion AeW^ + Awg = in O with the boundary condition S^w^ = on dO. 
Here S^Wg is the conormal derivative of w^, is some function defined in 
terms of the solution uq of the homogenized problem and containing rapidly 
oscillating coefficients (see §4). 

Some technical tricks, in particular, using the Steklov smoothing operator, 
are borrowed from [ZhPas]. 

First we prove estimate (0.7). For this, it is necessary to estimate the 
i?^-norm of by C£^/^||F||^2(0)- Next we obtain (0.6) by estimating the 
La-norm of in terms of C£||F 11^2(0)- For this, we rely on the (already 
proved) estimate (0.7) and use the duality arguments. 

0.6. Plan of the paper. The paper consists of nine sections. In §1, we 
introduce the class of operators in L2{W^; C") that we consider, describe the 
effective operator and the corrector, and give the results from [BSu2,4] and 
[PSu2] needed in what follows. Precisely, we deduce theorems about approx- 
imation of the resolvent {Ae + A/)~^ for arbitrary A > from the known 
results for the operator [Ac + I)~'. In §2, the statement of the Neumann 
problem in a bounded domain is discussed, some auxiliary material is given. 



HOMOGENIZATION OF THE NEUMANN PROBLEM FOR ELLIPTIC SYSTEMS 5 



and the homogenized problem is described. In §3, different auxihary state- 
ments are collected. Main results (Theorems 4.1 and 4.2) are formulated in 
§4. Also, in §4 the first part of the proof is given: the boundary layer correc- 
tion term is introduced, and the proofs of Theorems 4.1 and 4.2 are reduced 
to the estimates for this correction term. §5 contains the proof of Theorem 
4.2, and §6 contains the proof of Theorem 4.1. In §7 the case where A G -Lqo 
is studied; in this case it is possible to get rid of the smoothing operator 
in the corrector. In §8, estimates in a strictly interior subdomain O' of the 
domain O are obtained: it is shown that estimate (0.6) and the results in 
imply estimate (0.8). Finally, §9 is devoted to the case where A = 0: 
the corresponding results are deduced from the theorems for the case where 
A > 0. 

0.7. Notation. Let and ii* be complex separable Hilbert spaces. The 
symbols (•, •)/, and || • \\sj stand for the inner product and the norm in Sj; the 
symbol || • stands for the norm of a linear continuous operator acting 

from Sj to ij*. By / = If, we denote the identity operator in Sj. 

The symbols (•, •) and | • | stand for the inner product and the norm 
in C"; 1 = 1„ is the identity (n x n)-matrix. If a is an (n x n)-matrix, 
the symbol \a\ denotes the norm of a as an operator in C"". Wc denote 
X = {xi, . . . ,Xd) G K'^, iDj = dj = d/dxj, j = l,...,d, D = — iV = 
{Di, . . . ,D(i). The Lp-classes of C"-valued functions in a domain O C 
are denoted by Lp(0;C"), 1 ^ p ^ oo. The Sobolev classes of C"-valued 
functions in a domain O C M'' arc denoted by H'{0; C"). Next, H^{0; C") 
is the closure of C^{0; C") in H\0; C"). If n = 1, we write simply Lp{0), 
H^{0), etc., but sometimes we use such abbreviated notation also for the 
spaces of vector- valued or matrix- valued functions. 
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§1. HOMOGENIZATION PROBLEM FOR A PERIODIC ELLIPTIC OPERATOR IN 

L2(R'^;C") 

In this section, we describe the class of matrix elliptic operators that we 
consider and give the results for homogenization problem in obtained 
in [BSu2], [BSu4], and also in [PSu2]. Precisely, from the known results 
about approximation of the resolvent (As + we deduce theorems about 
approximation of the resolvent {Ae + XI)"^ for any A > 0. 
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1.1. Lattices in W^. Let F C M"^ be a lattice generated by the basis 
ai,...,ad gM'^': 

d 

F = {a G M'' : a = ^ i/^a^-, uj G Z}, 

i=i 

and let O be the (elementary) cell of the lattice F: 

d ^ ^ 

:= {x G R'^ : X = r,a,-, -- < < -}. 

We use the notation \^\ = measil. The basis hi, . . . ,hd in M.'^ dual to 
ai, . . . , a^ is defined by the relations (bj, a^) = 2iT5ij. This basis generates 
a lattice F dual to the lattice F. Below we use the notation 

ro = - min_ |b|, ri = -diamri. 
2 o^ber 2 

By we denote the subspace of all functions in H^{Q,) whose F- 

periodic extension to belongs to H^^^CR'^). If (^(x) is a F-periodic function 
in W^, we denote 

(p^{x) := v?(e"^x), e > 0. 

1.2. The class of operators. In L2(M°';C"), consider a second order DO 
Ae formally given by the differential expression 

A = 6(D)*/(x)6(D), e>0. (1.1) 

Here ^(x) is a measurable Hermitian (m x m)-matrix-valued function (in 
general, with complex entries). We assume that gix) is periodic with respect 
to the lattice F, bounded and uniformly positive definite, i. e., 

dm ^ 5(x) ^ cl m, xgM; 0<c^c< go. 

Next, b(D) is a homogeneous (m x n)-matrix first order DO with constant 
coefficients: 

d 

b{Ty) = J2blDi. (1.2) 

1=1 

Here bi are constant matrices (in general, with complex entries). The symbol 
b{^) = Yl'i=iblCh ^ £ I^'^j corresponds to the operator 6(D). It is assumed 
that m'^ n and 

rank b{^) = n, Q^^eR'^. (1.3) 
This condition is equivalent to the following inequalities 

aoln ^ b{e)*b{e) ^ ailn, 9 G S'^'S < ao ^ ai < oo, (1.4) 
with some positive constants ao and ai. Prom (1.4) it follows that 

\h\^a]''', l = l,...,d. (1.5) 
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The precise definition of the operator Ag is given in terms of the corre- 
sponding quadratic form 

ae[u,u\= [ (/(x)6(D)u,6(D)u) dx, u e H\R'^;C'). 

Under the above assumptions this form is closed in L2(M'^; C") and nonneg- 
ative. Using the Fourier transformation and condition (1.4), it is easy to 
check that 

Co / iDupdx ^ a£[u,u] ^ ci /" |Dupdx, u G i?^(M'^; C"), (1.6) 

where cq = aolb^'^IIZi' '^i = "ilblU^o- 

The simplest example of the operator (1.1) is the scalar elliptic operator 
Ae = — div3^(x)V = D*g''^(x)D. In this case we have n = 1, m = d, 
6(D) = D. Obviously, condition (1.4) is valid with = ai = 1. Another 
example is the operator of elasticity theory which can be written in the 
form (1.1) with n = d and m = d{d+ l)/2. These and other examples are 
discussed in [BSu2] in detail. 

1.3. The effective operator. In order to formulate the results, we need 
to introduce the effective operator A^. 

Let an (n x m)-matrix-valued function A(x) be the (weak) F-periodic 
solution of the problem 

KD)*5(x) (KD)A(x) + 1„0 = 0, / A(x) dx = 0. (1.7) 

Jn 

In other words, for the columns Vj (x), j = 1, . . . , m, of the matrix A(x) the 
following is true: Vj G H^{Q,;C"'), we have 

/ (5(x)(6(D)v,(x)+%),6(D)r7(x))dx = 0, Vr? G ^^(J^; C"), 
Jn 

and Vj(x) dx = 0. Here ei, . . . , e^^ is the standard orthonormal basis in 
C™. 

The so-called effective matrix of size m x m is defined as follows: 

g'' = \n\-' [ 9(x)dx, (1.8) 
Jn 

where 

^(x):=5(x)(6(D)A(x) + l^). (1.9) 

It turns out that the matrix is positive definite. The effective operator 
A^ for the operator (1.1) is given by the differential expression 

A^ = 6(D)*5°6(D) 

on the domain if2(M"f; C"). 

Below we need the following estimates for A(x) proved in [BSu3, (6.28) 
and Subsection 7.3]: 

||DA|U,(n) ^ M'/'m'/W%\\]!l\\9-rt (1-10) 
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1.4. Properties of the effective matrix. The following properties of the 
effective matrix are proved in [BSu2, Chapter 3, Theorem 1.5]. 

Proposition 1.1. The effective matrix satisfies the estimates 

5 ^ 5° ^ 5- (1-12) 

Here 

If m = n, the effective matrix coincides with g. 

In homogenization theory for specific DO's, estimates (1.12) are known 
as the Voight-Reuss bracketing. We distinguish the cases where one of the 
inequalities in (1.12) becomes an identity. The following statements were 
obtained in [BSu2, Chapter 3, Propositions 1.6 and 1.7]. 

Proposition 1.2. The identity g^ = g is equivalent to the relations 

6(D)*gfc(x) = 0, k = l,...,m, (1.13) 

where gfc(x), k = 1, . . . ,m, are the columns of the matrix 5'(x). 
Proposition 1.3. The identity g^ = g is equivalent to the relations 

lfe(x) = lO + 6(D)wfe, l^GC"*, WkeH\n-Cn, k = l,...,m, (1.14) 
where Ijk(x), k = 1, . . . ,m, are the columns of the matrix g{y:)~^. 

Obviously, (1.12) implies the following estimates for the norms of the 
matrices g^ and {g^)~^'- 

\9'\^\\9\\l^, \{9'r'\^\\9-'\\L^- (1-15) 

Note that, by (1.4) and (1.15), the symbol of the effective operator 
satisfies the following inequality: 

6(0 VKO > co|^|2l„, ^GM'^, CO = ao\\g-^\\^\ (1.16) 

1.5. Smoothing in Steklov's sense. Let be the operator in 

L2(M'^;C'") given by 

{SeU){x) = [ U(x - £Z) dz (1.17) 

and called the Steklov smoothing operator. Note that 

\\^e\\L2(M'i)^L2(m'^) ^ 1- (1-18) 

Obviously, D°'SeU = SsD^'n for u G H'iR'^; C™) and \a\ ^ s. Therefore, 

ll'S'e||/fs(iRd)^_f/s(I[jd) ^1, s G N. (1-19) 



We mention some properties of the operator (1.17), see [ZhPas, Lemmas 
1.1 and 1.2] or [PSu2, Propositions 3.1, 3.2]. 
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Proposition 1.4. For any u G H\R'';C"') we have 

ll'S'eU - u||i2(K<*) ^ e''l||Du||i2(Kd), £ > 0. 

Proposition 1.5. Let /(x) be a T-periodic function in such that f G 
L2{'^). Let [/^] be the operator of multiplication by the function /'^(x). Then 
the operator [f^jSi. is continuous in L2(M'*;C™), and 

\\[f]Se\\L2{Rd)^L2(R'i} ^ l^r"^^^ll/l|L2(n), £ > 0. 

1.6. The results for homogenization problem in W^. Consider the 
following elliptic equation in W^: 

AeUe + AUe = F, (1.20) 

where F G L2(K'^;C"^), and A > is a parameter. It is known that the 
solution Ue converges in L2(M'^;C") to the solution of the "homogenized" 
equation 

^Ouo + Auo = F, (1.21) 

as £ ^ 0. 

Theorem 1.6. Let Ug be the solution of the equation (1.20), and let Uq be 
the solution of the equation (1.21). Then 

l|U£ ~ Uo||l2(K=';C") ^ C'i(A)£||F||^2(R=';C")) £ > 0, 

or, in operator terms, 

||(.4e + A/)-i-(.4° + A7)-i||i^(Kd;cn)_,L,(Mrf;Cn) ^ C'i(A)£, £>0. (1.22) 

Here Ci(A) = CiX^^^^, and the constant Ci depends only on the norms 
IblUoo) lb~^llLoo5 constants ao, ai from (1.4) and the parameters of 

the lattice T. 

Proof. In [BSu2, Chapter 4, Theorem 2.1], estimate (1.22) was proved in 
the case where A = 1 for < e ^ 1. We only have to explain how this 
estimate is carried over to the general case. 

Note that for A = 1 and e > 1 the left-hand side of (1.22) is obviously less 
than 2, and then also less than 2e. Therefore, wc start with the estimate 

II (^, + I)-i - (^0 + /)-iL,(ffid)^L,(Md) ^ C,e, s > 0. (1.23) 

Here the constant Ci depends only on H^Hloo, ||5~''"||loo> ^0, (Xi, and the 
parameters of the lattice T. 

Next, by the scaling transformation, (1.23) is equivalent to the inequality 

\\(A + e^I)-'-{A'' + e^I)-'h^^^,^^L^^^,)^ Cie-\ e > 0. (1.24) 

Here A = b{D)*g{x.)b{D). Using (1.24) with e replaced by eA^/^ and ap- 
plying the inverse transformation, we arrive at (1.22) with the constant 
Ci(A) = (7iA-V2. . 
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In order to approximate Ug in H^{W^; C") it is necessary to take the first 
order corrector into account. We put 

K^{s) = [A^]5,6(D)(^o + Xir\ (1.25) 

Here [A^] is the operator of multipUcation by the matrix-valued function 

A(e~^x), and is the smoothing operator defined by (1.17). The operator 
(1.25) is a continuous mapping of L2(M'^;C") into H\R'^;C"-). Indeed, the 
operator b{'D){A^+\I)-^ is continuous from L2(M'^; C") to H^R'^; C™). The 
operator [A^]Se is a continuous mapping of H^{R'^; V^) into H^iW^; C"). It 
can be easily checked by using Proposition 1.5 and relation A G H^[Q). Wc 
have e\\Kx{e)\\]^^_^fji = 0(1) for small e. (See the proof of Theorem 1.7 
below, where this is checked for A = 1.) 

"The first order approximation" to the solution is given by 

= uo + eA^S,b(D)uo = {A^ + Xiy^F + eKx{e)F. (1.26) 

Theorem 1.7. Let Ug be the solution of the equation (1.20), and let uq 
be the solution of the equation (1.21). Suppose that is defined by (1.26). 
Then 

ll^e ~ Vell^i^jjd.Cn) ^ C'2(A)£||F||^2(]Rd;C")> £ > 0, (1-27) 

or, in operator terms, 

||(.A, + A/)-^-(.4° + A/)-^-£KA(£)||L2(Md;C"Hffi(ffid;C") ^ C'2(A)e, e>0. 

Here C2(A) = C2(A~^/^ + 1), and the constant C2 depends only on 
d,m,ao,ai, \\g\\Looi ll5~^IUoo' "^^^ parameters of the lattice F. 
Proof. In [BSu4, Theorem 10.6], a similar result was proved for A = 1, but 
with a different smoothing operator instead of S^- In [PSu2, Theorem 3.3], 
it was shown that it is possible to pass to the smoothing operator S^, and 

(1.28) was proved for A = 1 and < £ ^ 1. We only have to explain how 
this estimate can be carried over to the general case. 

Thus, we start with the estimate 

\\iAe + ir'-{A^ + ir'-eK,{e)\\L,(^^a)^HHRi)^ C'2£, < £ ^ 1. (1.29) 

The constant C2 depends only on d, m, ag, ai, II^IIl^^ , \\g~^ \\loo > the pa- 
rameters of the lattice F. 

For £ > 1 estimates are trivial: each operator under the norm sign in 

(1.29) is estimated separately. Prom the lower inequality (1.6) it follows 
that 

min{co,l}||(.4e + /)"V||^i(^d) ^ {{^6 + I)~'^y,y)L2{Rd) ^ Ml^(Rd) 

for V G L2(R'';C"), whence 

WiAe + ly^h^iRd^^miR") ^ max{l,Co^/^} = max{l,ao^''^||5"i||^/^}. 

(1.30) 
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By (1.15), the norm of {A^ + I)^^ satisfies the same estimate: 

WiA^ + I)''^\\L2(R'i)^m{R'i) ^ max{l,Co^/^}. (1.31) 

Now we estimate the (L2 — )■ H^)-noicm of the operator £Ki{e) = 
e[A^]Seb{'D){A^ + I)-^. Let F e L2(M'^;C"). Then, by Proposition 1.5, 

||eKi(e)F||i^(K.) ^e|J^|-^/2||A|U,(f,)||6(D)(^o + /)-iF||^^(K.). 
Using the Fourier transformation and (1.4), (1.16), we obtain: 

\\b{B){A' + I)-'F\\l^^,^ ^ [ |6(^)|2|(6(^)V6(0 + l)-Y|F(Opd^ 
^ ai |^P(co|^P + l)-'m)\' ^ a,{2co)-'\\F\\l^^,y 

(1.32) 

Here F(^) is the Fourier image of F(x). Hence, 

||eKi(e)F||i^(K.) ^ e|Q|-^/2||^||^^^^^^i/2^2co)-i/2||p||^^^^_^^_ (1_33^ 
Consider the derivatives 

edj(Ki{e)F) = [{djAY]Seb{B){A^ + I)-^F + e[A']Seb{B)dj{A° + ly^F. 
By Proposition 1.5, we have: 

J^||£5,-(Ki(£)F)||i^(^,) ^ 2|fi|-i||DA||i^(^)||6(D)(^o + J)-iF 



|2 



(1.34) 



Similarly to (1.32), 

Y,\m)dj{A' + ir'n 

j=i (1.35) 
^ ai^^ |^|'(co|^|' + l)-'m)fd^ ^ «iCo^||F||2^(^,). 

As a result, from (1.32), (1.34), and (1.35) it fohows that 
d 

^|0|-i||DA|||^(f,)aiCoiF|| 



-2ir>l-l|l A Il2 „ „-2in7i|2 



(1.36) 



+ 2£^|J^|-lA||i^(^)aiCo-^||F||i^(^,). 
By (1.33) and (1.36), 

\\eKi{e)\\L,(j^a)^miR'i) ^ M'^'' al^'^c^' /''\\BAU^^n) 

+ e|OrV2«f||A||,^(^^(2^)-V2(i + 4,-1)1/2. 
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Combining this with (1.10) and (1.11), we obtain 

lk^l(^)llL2(K'^)^J?HIK'*) ^^3 + C4£, £ > 0, (1.37) 

where 

^ 1/2 1/2 -111 ||1/2|| -111 

Relations (1.30), (1.31), and (1.37) imply that 

\\{Ae + I)-' ~ (^° + /)-' -£i^l(e)||L,(R^H//i(R<*) 

^ 2max{l,ao^/^||5-^||};'^} + C3 + C4£, £ > 0. 

Obviously, for e > 1 the right-hand side does not exceed C2£, where C2 = 
2max{l,Q!Q ^^^Ib^^ll^^f } + C3 + C4. Combining this with (1.29), we obtain 
estimate of the form (1.29) for all £ > 0: 

IKA + /)"'-(^° +/)-'- £i^l(£)||L2(RdHi^i(M<^) ^C*2£, £>0, (1.38) 

where C2 = max{C'2, C2}. 

A similar approximation for the operator (Ae+XI)'^ with arbitrary A > 
can be easily deduced from (1.38). Indeed, (1.38) is equivalent to 

||(_A + /)V2 ((^^ + J)-l _ (^0 + ^)-l _ ^^^(^)) ||^,(K.)^L,(M.) C2£ 

(1.39) 

for £ > 0. By the scaling transformation, (1.39) is equivalent to 
||(-A + e2/)V2 ((^ + £2^)-i _ (^0 +£2/)-! _ A56(D)(^o + £2/)-i) \\^^^^^ 

^C2, £ > 0. 

(1.40) 

Here 5 = 5i is the operator (1.17) with £ = 1. From (1.40) with e replaced 
by £A^/^, applying the inverse transformation, we obtain: 

||(-A + A/)V2 ((^^ + A/)-i - (^0 + A7)-i - eK^ie)) ||i,(M.)^L,(M.) ^ ^2£ 

for £ > 0. This implies (1.28) with C2(A) = CiiX'^^^ + 1). • 

Theorem 1.7 allows one to obtain approximation for the flux := 
/6(D)ue in L2(M'^;C"^). 

Theorem 1.8. Let be the solution of the equation (1.20), and let uq be 
the solution of the equation (1.21). Let Pe := (/^6(D)u£. Then we have 

\\Pe - ff^S'e6(D)uo||i2(i;d.cm) ^ C5(A)£||F||i2(i:'';C")> £ > 0- (1-41) 

Here g{x) is the matrix (1.9), C^iX) = C'^X~^^'^ + ^'5') ^'^'^ the constants C'^, 
C5 depend only on d, m, ao,ai, II^Hloo' lb~"^lkoo' parameters of the 

lattice r. 

Proof. By (1.4) and (1.27), 

Hp, - /6(D)v,||i,(i,d) ^ £a}/'||5||L^C2(A)||F||^^(Kd), e > 0. (1.42) 
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From (1.2) and (1.26) it follows that 

d 

/6(D)v, = /6(D)uo + /(6(D)A)^S,6(D)uo + £ J]/6jA^5,6(D)Auo. 

1=1 

(1.43) 

Using (1.5) and Proposition 1.5, we estimate the last term in the right-hand 
side of (1.43): 

d 

^g%A'S,bi'D)Dino 

1=1 (1-44) 

^e||5l|L^ai^'l^^r'/'l|A||L2(Q)d'/'l|D6(D)uo||i,(K.). 
Next, by Proposition 1.4, 

||/6(D)uo - /S'£6(D)uo||L2(Rd) ^ e||5||L^ri||D6(D)uo||i2(Kd)- (1-45) 
By (1.9), we have 

/5,6(D)uo + /(6(D)A)^S,6(D)uo = rSeb(D)uo. (1.46) 
Combining (1.43)-(1.46) and (1.11), wc arrive at 

||/6(D)v,-5^5,6(D)uo||i,(Md) <C6e||D6(D)uo||i,(Kd), e > 0, (1.47) 

where = (dm)^/2(2ro)-ia}''^Q;Q ^ri\\g\\L^. 
Similarly to (1.35), 

||D6(D)uo||i^(K.) ^ m j^^ \i\\c^\$,? + A)-2|F(4)P di ^ a^cfml^^^ay 

(1.48) 

Combining (1.42) and (1.47), (1.48), we obtain (1.41) with the constant 
C5(A) = aA-V2 + ^i^g^g ^ C2a]'^\\g\\L^, C'i = C2ay^\\g\\L^ + 

^ 1/2 -111 -111 

Now we distinguish the case where the corrector is equal to zero. Next 
statement follows from Theorem 1.7, Proposition 1.2, and equation (1.7). 

Proposition 1.9. Let Ug he the solution of the equation (1.20), and let Uq 
he the solution of the equation (1.21). If g^ =9, i- e., relations (1.13) are 
satisfied, then A = 0, Kx{e) = 0, and we have 

ll^e - uo||//i(Rd;C") ^ C2{\)e\\F\\i^^(^d), e > 0. 



1.7. The results for homogenization problem in R'^ in the case 
where A G L^o- It turns out that under some additional assumptions on 
the solution of the problem (1.7) the smoothing operator in (1.25) can 
be removed (replaced by the identity operator). We impose the following 
condition. 

Condition 1.10. Suppose that the T-periodic solution A(x) of the prohlem 
(1.7) is bounded: A G Loo- 
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We need the following multiplicative property of A, see [PSu2, Corollary 
2.4]. 

Proposition 1.11. Under Condition 1.10 /or any functionu G H^{M.^) and 
e > we have 

[ \{J^A)%:K)f\u\U:K^l3i\\u\\l.^a)+P2\\M\iyf |D«pdx. 

The constants Pi, j32 are given by 

/3i = 16maQ^\\g\\L^\\g-^\\L^, 

/32 = 2(1 + 2dao ^ai + 20dao'ai\\g\\L^ \\g-^ 

We put 

i^^(£) = [A^]6(D)(^o + A/)-i. 

By Proposition 1.11, it is easily seen that under Condition 1.10 the operator 
K'^ie) is a continuous mapping of L2(M'^;C") to ifi(R'^;C"). Instead of 
(1.26), consider another first order approximation of u^: 

= uo + eA^5(D)uo = (^° + Xiy^F + eK'i{e)F. (1.49) 

Theorem 1.12. Suppose that Condition 1.10 is satisfied. Let be the 
solution of the equation (1.20), and let Uq be the solution of the equation 
(1.21). Let Ve be defined by (1.49). Then 

W^e - ^e\\m(Rd-C") ^ C'7(A)e||F||^2(K'i;C")' £ > 0, (1.50) 
or, in operator terms, 

||(^, + A/)-i-(^o + A/)-^-£K°(£)||i^(Kd;CnH^^i(ffid;C") ^ CrWs, e > 0. 
For the flux Pe = g^b(D)us we have 

\\Pe - 5^b(D)uo||i2(Rd;C-) ^ C'8(A)£||F (IR'';C») > ^ > 0> (1-51) 

where g{x) is the matrix (1.9), Cr{\) = C'^X'^I'^ + C^', C^{\) = 
CgA"^/^ + Cg, and the constants C^, C", C'g, Cg depend only on 
d,m,ao,ai, \\g\\Loo^ Ib^^lUooJ parameters of the lattice F, and the norm 

I|A||loo- 

Proof. For the proof of (1.50) we need to estimate the iJ^-norm of the 
function eA^{L — Ss)b{'D)uo- We start with the L2-norm. By Condition 1.10 
and (1.18), 

£||A^(7-5,)6(D)uo||i,(Md) ^ 2£||A||i^||6(D)uo||i,(Kd). (1-52) 

Consider the derivatives 

di{eA'{I - 5,)6(D)uo) = (5;A)^(/ - 5,)6(D)uo + eA'{L - 5,)6(D)5;Uo. 

(1.53) 
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The second term in the right-hand side of (1.53) is estimated with the help 
of Condition 1.10 and (1.18): 

d 

j;i|eA-(/-5,)6(D)aiUo||i^(^.) ^4£2||A||2^||D6(D)uo||i^(jj.). (1.54) 

(=1 

Using Proposition 1.11, we estimate the first term in the right-hand side of 
(1.53): 

d 

^||(W(/-S,)6(D)uo||i^(i,.) 

1=1 

d 

^ ^i||(/- 5,)6(D)uo||i^(K.) +/32£'||A||L E - Se)h{Ti)Din4l^^^,y 

1=1 

(1.55) 

By Proposition 1.4, 

11(7 - 5e)6(D)uoL,(ffid) ^ £ri||D6(D)uo||. (1.56) 

Now from (1.55), (1.56), and (1.18) it follows that 

d 

1=1 

(1.57) 

As a result, combining (1.53), (1.54), and (1.57), we arrive at 

d 

Y,\m^'{i-Se)hm "0)|li2(Rd) 

1=1 (,i.58j 
^ (4(1 + ;32)||A||L +;3ir?) ||D6(D)uo||i^(Kd). 
Similarly to (1.32), we have 

||6(D)uo||i^(j,,) ^ ai |^|2(co|^|2+A)-2|F(0|'d^ ^ ai(2Aco)-iF||2^( 

Together with (1.52), (1.58), and (1.48) this yields 

||£A^(J - S,)6(D)uo||^i(Kd) ^ C9(A)£||F||^2(Kd), (1.59) 
where C9(A) = C^A-V2 + ^ C, = V2||A|U^a}/'ao 

= ^/2af ao ^ Wg-' (2(1 + /32)^/^ || A||^^ + /3j/Vi) . 



As a result, the desired estimate (1.50) follows from (1.26), (1-27), (1.49), 
and (1.59); the constant CriX) is given by CriX) = C^A'^/^ + where 

C'^ = C2 + Cg, = C2 + Cg . 

To prove (1.51), note that from (1.50) and (1.4) it follows that 

||p,-/6(D)v,||^,(K<i) ^£C7(A)||5||L^«}/'||F||i2(K<i). (1.60) 
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Next, by (1.2), (1.9), and (1.49) 

d 

/6(D)ve = 5'&(D)uo + e ^ g%A%{B)DiUo. (1.61) 

1=1 

The last term in the right-hand side of (1.61) is estimated with the help of 
(1.5), (1.48), and Condition 1.10: 

d 



(=1 



^ ^IblUoollAIU^aiCo V'/'||F||i^(K.). (1.62) 

Now relations (1.60)-(1.62) imply (1.51) with the constant Cs{X) = 
C^A-V2 + C'l where C, = O^WgU^al^^ 

In some cases Condition 1.10 is valid automatically. Next statement was 
proved in [BSu4, Lemma 8.7]. 

Proposition 1.13. Condition 1.10 is a fortiori valid if at least one of the 
following assumptions is fulfilled: 
r. d ^ 2; 

2°. d ^ 1 and Ae = D*gi^(x)D, where g{x) has real entries; 

3°. dimension is arbitrary, and g^ = g, i. e., relations (1.14) are satisfied. 

Note that Condition 1.10 is also ensured if ^(x) is smooth enough. 

We distinguish the special case where g^ = g- In this case the matrix 
(1.9) is constant: (/(x) = g^ = g; moreover, Condition 1.10 is satisfied. 
Applying the statement of Theorem 1.12 concerning the fluxes, we arrive at 
the following statement. 

Proposition 1.14. Suppose that is the solution of the equation (1.20), 
uo is the solution of the equation (1-21), and p^ = 5^6(D)u£. Let g^ = g, 
i. e., relations (1.14) are satisfied. Then we have 

IIPe - 5°^(D)uo||i2(Kd.cm) ^ C'8(A)£||F||i2(Md;C")> £ > 0. 

§2. The Neumann problem in a bounded domain: preliminaries 

2.1. Coercivity. Let O C M'' be a bounded domain with the boundary 
of class C^'^. We impose an additional condition on the symbol h{$,) = 

Ef=i h^l for C G C^. 

Condition 2.1. The matrix-valued function b{^), ^ G C^, is such that 

rank6(0 = n, 0^ie&. (2.1) 

Note that Condition 2.1 is more restrictive than (1.3). According to 
[Ne] (see Theorem 7.8 in §3.7), Condition 2.1 is necessary and sufficient for 
coercivity of the form ||5(D)u|||^^q^ on H^{0; C") (moreover, this is true in 
any bounded domain O with the Lipschitz boundary). 
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Proposition 2.2. [Ne] Condition 2.1 is necessary and sufficient for exis- 
tence of constants Ci > and C2 ^ such that the Carding type inequality 

IIKD)"llL(0)+C2||u||i^(^) >Ci||Du||i^(0), ueH\0;C^), (2.2) 

is satisfied. 

Remark 2.3. The constants Ci and C2 depend on the matrix b{^) and the 
domain O, but in general case it is difficult to control these constants ex- 
plicitly. However, often for particular operators they are known. Therefore, 
in what follows we indicate the dependence of other constants on Ci and C2. 

In what follows, we assume that Condition 2.1 is satisfied. 
2.2. Statement of the problem. In L2(0;C"^), consider the operator 
AN,e formally given by the differential expression 6(D)*g'^(x)6(D) with the 
Neumann condition on dO. The precise definition of An,s is given in terms 
of the quadratic form 

a^,,[u,u] := / (/(x)6(D)u,6(D)u) dx, ueH\0;C''). (2.3) 
Jo 

By (1.2) and (1.5), 

ajv,.[u,u] ^dai||5||L^||Du||i^(0), u G /^^(O; C"). (2.4) 
Prom (2.2) it follows that 

a^,,[u,u] > ll^-^Z^ (Ci||Du|||^(o)-C2||u||i^(0)) , ueHHO;Cn. 

(2.5) 

By (2.4) and (2.5), the form (2.3) is closed and lower semibounded. By 
definition, AN,e is the selfadjoint operator in L2(0;C^) generated by the 
quadratic form (2.3). 

Let X> be a parameter subject to the following restriction: 

>^>C2\\9-'\\ll. (2-6) 

Then (2.5) implies that 

a7v,e[u,u] + A||u|||,^(0) ^ ca||u|||i(0), uGii"\0;C"), 

cx :=mm{C,\\g-Xl,\-C2\\g-%l}. ^^'^ 

Hence, the operator An,£ + A/ is positive definite. 

Our goal is to approximate the inverse operator {AN,e + A/)~^ for small 
£ in the operator norm in L2(0;C") and in the norm of operators acting 
from L2(0;C") to ii"^(0;C"). In terms of solutions, we are interested in 
the behavior of the generalized solution G H^{0;C^) of the Neumann 
problem 

6(D)*/(x)6(D)u,(x) + Au,(x) = F(x), x G O; 9>,|50 = 0, (2.8) 

where F G ^2(0; C"). Then = {AN,e + A/)-^F. 

Here the notation for the "conormal derivative" was used. Let i/(x) = 
Yl'i=i be the unit outer normal vector to dO at the point x G dO. 
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Here ei,...,ed is the standard orthonormal basis in W^. Introduce the 
matrix 6*(i/(x)) = Yli=i^^l{'^)- Then formally the conormal derivative is 
given by 

d 

9tu(x) := 6*(iy(x))/(x)6(9)u(x) = ^ b^g'{^)bM^)9M^)^ ^ ^ 90. 

Note that, under our assumptions, (2.8) has formal sense. By definition, 
we say that € H^{0;C^) is the generalized solution of the problem (2.8) 
if Ue satisfies the following integral identity 

I ((/6(D)u„6(D)r7) + A(u„77)) cix= /(F,77)dx, ^r^ e H\0-C^). 
Jo Jo 

(2.9) 

Taking (2.4) and (2.7) into account, we can view the form ajv^£[u,T7] + 
A(u, J7)£,2(C') as a (new) inner product in i^^(0;C"). The right-hand side of 
(2.9) is an antilinear continuous functional of G i7^(0;C"). By the Riss 
theorem, there exists unique solution u^; by (2.7), it satisfies the estimate 

lluelkHO) ^ c-^I|F||l2(0)- (2-10) 
In operator terms, (2.10) means that 

\\{AN,e + >^lT^\\L^{^0)^H^{0) < C-\ (2.11) 

2.3. The trace and extension operators. It is well known that, under 
our assumption on the domain O (that the boundary is of class C^'^), the 
trace operator 7 taking a function u in O into its trace on the boundary dO 
is defined correctly linear continuous operator 

7 : H'{0- C") ^ H'-^l^{0- C"), s = 1, 2. (2.12) 

Herewith, 

^^1' ^2 13) 

H'^(0)^H^/'^{dO) ^ "^2, 

where the constants ci and C2 depend only on the domain O. 

There exists a linear continuous extension operator taking a function at 
the boundary into its extension to the domain O (this is the right inverse 
to the operator (2.12)). Obviously, such an operator is not unique. It 
is convenient to choose this operator in the following way. Consider the 
Dirichlet problem 

-Au-|-u = in O, ulao = (2.14) 

where if G H^/'^{dO; C"^). It is well known that there exists a unique solution 
u G H^{0; C"), and it satisfies the estimate ||u||j:i-i(0) ^ ci||<p||j|^i/2(gQ^. The 
constant ci depends only on the domain O. Denote by T the operator taking 
a function (p into the solution u of the problem (2.14): u = Tcp. Then 

T : H^/'^{dO; C") ^ H^{0; C") (2.15) 
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is a linear continuous operator, T is the right inverse of the operator (2.12) 
(with s = 1), and 

ll^lli?i/2(ao)^j?i(0) ^ ci. (2-16) 

By the theorem about regularity of solutions (see, e. g., [McL, Chapter 
4]), if ipe H^/^{dO-C''), then u G F2(C);C") and the operator 

T : H^/^{dO; C") ^ H^{0; C") (2.17) 

is continuous. Its norm satisfies the estimate 

1 1 '^11^3/2 (90)^^2(0) ^ C2, (2.18) 

where the constant C2 depends only on the domain O. The operator (2.17) 
is the right inverse of the operator (2.12) with s = 2. 

2.4. Definition of the conormal derivative. The Neumann problem 

for the homogeneous equation. We need to give the precise definition of 

the conormal derivative (corresponding to the operator As). The conormal 
derivative df,{ is defined as an element of H~^^^(dO; C"), definition is based 
on the Green formula, in which the function f G i7^(0;C") and the right- 
hand side of the equation arc given. For our purposes, it suffices to give the 
definition of the conormal derivative for a function f G i?^(C;C") in the 
case where Agf belongs to L2{0; C"). 

Thus, let f G H^iO:C'') and G L2(C;C") be such that Aef = 
inside O, which is understood as the integral identity 

/(/6(D)f,6(D)77)dx= /(*„r7)dx, ^rj e H^{0;C^). (2.19) 
Jo JO 

Let le[f] be the following antilinear functional of (p €z H^/'^{dO;C^): 

leW>\= I (/6(D)f,6(D)u)c?x- / (*e,u)c?x, (2.20) 
Jo JO 

where u G H^{0;C^) is a function such that 7U = (p. By (2.19), the right- 
hand side of (2.20) does not depend on the choice of an extension u of the 
function ip\ therefore, definition of the functional (2.20) is correct. It is 
convenient to put u = Tcp, where the operator T is defined in Subsection 
2.3. Then, by (1.2), (1.5), and (2.16), the functional (2.20) is continuous: 

|ZeMI ^ ||*£||l2(0)||u||l2(o) + aid||5||L^||Df||i2(o)||Du||i2(0) 21) 

< Ci {\\^e\\L.2{0) + (^ld\\9\\Loo¥\\m{0)) 1 1 ¥5 1 1 if 1/2 (ao)- 

Recall that the space iJ-V2(ac'; c") is dual to ijV2(5e). c") ^jth re- 
spect to the L2(9C'; C")-duality. It is common to write the meaning of 
the functional ^ € H-^l'^{dO;([7') on the element G H^/'^{dO;C'^) as 
(-j/j. If) L2{dO) (which is extension of the inner product in L2{dO; C") for the 
pairs in iJ-V2(ac'; C") x H^/'^{dO-C^)). Herewith, 

l(V>,</g)L.(a0)l , . 

H-y^idO) - sup -J— i . (2.22) 

o^<^ei?i/2(ao) 1 1 VI I if 1/2 (ao) 
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For the antilinear continuous functional le['^\ on H^^^{dO; C") defined by 
(2.20) there exists a unique element e H~^/^{dO;C'") such that 

k[^] = {iP„^)L^^90), V>&H^^\dO;Cn. (2.23) 

By definition, we say that if)^ is the conormal derivative of a function f and 
write = df,f. From (2.21)-(2.23) it follows that 

||5^f 11^-1/2(90) ^ Ci (W^ehiiO) + "I'^lblUoc l|f • 

Now we discuss the statement of the Neumann problem for the ho- 
mogeneous equation with the nonhomogeneous boundary condition. Let 
■0 G H~^/'^{dO;C^). Consider the generalized solution of the problem 

Aere + Xre = (i in O; dlr^ = ip on 90. (2.24) 

By definition, the generalized solution of the problem (2.24) is an element 
G H^{0;C"') satisfying the identity 

((/6(D)r„6(D)r7) + \{v,,r^)) dx = (V',7r7)L.(ao), Vt? G H\0-C'). 

(2.25) 

The following statement is checked in a standard way. 

Proposition 2.4. Let ijj G if^^/^(50; C"). Then the generalized solution 
G H^{0\C^) of the problem (2.24) exists, it is unique and satisfies the 
estimate 

^ c^^ci||i/>||^-i/2(9(p.cn), (2.26) 
where the constant c\ is defined by (2.7), andci by (2.13). 

Proof. As has been already mentioned, the form in the left-hand side of 
(2.25) can be viewed as an inner product in i7-'^(C'; C"). By (2.13) and 
(2.22), the right-hand side of (2.25) is an antilinear continuous functional of 
77 G i^HC; C"), and we have 

I('0,7^7)l2(O)I «S ci||'0||j:^-i/2(aci)||r7||j/i(ci), rj £ H^{0;C). 

Then, by the Riss theorem, the solution exists and is unique. By (2.7), 
it satisfies the estimate (2.26). • 

2.5. The "homogenized" problem. In L2(0;C"), consider the selfad- 
joint operator A% generated by the quadratic form 

a%[u,u]= [ {g%(D)u,b(D)u) dx, ueH\0;C). (2.27) 
Jo 

Here is the effective matrix defined by (1.8). By (1.15), estimates (2.4), 
(2.5) remain true for the form (2.27). 

Suppose that A still satisfies (2.6). Let uq G H^{0; C") be the generalized 
solution of the Neumann problem 

6(D)V6(D)uo(x) + Auo(x) = F(x), x G O; a°uo|ao = 0, (2.28) 
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where F G L2(C';C"). Then uq = {A% + A/)-^F. Here the formal differen- 
tial expression for the conormal derivative (corresponding to .4°) is of the 
form 

d 

5°u(x) = 6*(i/(x))505(9)u(x) = ^ b^g%,yii^)dM^), x G dO. (2.29) 

We say that Uq G i?^(0;C") is the generalized solution of the problem 
(2.28) if it satisfies the integral identity 

[ ((506(D)uo,6(D)r7) + A(uo,r?)) dx= /(F,77)dx, yrj G H\0;C''). 
Jo Jo 

(2.30) 

The solution exists, it is unique and satisfies the estimate 

l|uo||i/i(o) ^ c^^I|F||l2(0)- (2-31) 
In operator terms, (2.31) means that 

II (^^ + \I)-'\\l,(0)^hHO) ^ (2.32) 

By the condition dO G C^'^, the solution uq of the problem (2.28) satisfies 
uo G i?2(C);C") and 

l|uo||i/2(c).cn) ^ CA||F||i2(C>;C")- (2.33) 

Here the constant C\ depends only on the constants Ci and C2 from (2.2), on 
aO) Oil, IIs'IIloo' ll5~^llioo) ^1 the domain O. To justify this fact, we note 
that the operator hiY))* g^hiY)) is a strongly elliptic matrix operator with 
constant coefficients and refer to the theorems about regularity of solutions 
for strongly elliptic systems (see, e. g., [McL, Chapter 4]). 

From what was said it follows that the left-hand side of (2.28) belongs 
to L2(0;C"); equation is valid almost everywhere in O. The conormal 
derivative d^UQ is defined correctly by (2.29) as an element of i/^/^^^O; C"); 
the boundary condition in (2.28) can be understood in the sense of the trace 
theorem. Estimate (2.33) can be written as 

\\{A% + >^I)~^\\l2{0;C")^H^{0;€.-^) ^ C\. (2.34) 

In what follows, it is shown that the solution Ug of the problem (2.8) 
converges in L2{0\ C") to the solution uq of the homogenized problem (2.28), 
as £ ^ 0. We will estimate the error term ||u£ — Uq 11^2(0)- Also, we will find 
approximation of in H^{0\'C^). 

§3. Auxiliary statements 
This section contains various auxiliary statements needed below. 

3.1. Next statement is a version of the traditional lemma used in homoge- 
nization theory (see, e. g., [ZhKO, Chapter 1, §1]). 
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Lemma 3.1. Let fi{x), I = 1, . . . ,d, be T-periodic (n x m) -matrix-valued 
functions in such that 

„ d 

fieL2in), / /Kx)dx = 0, l = l,...,d; Va//(x) = 0, (3.1) 

where the last equation is understood in the distribution sense. There exist 
T-periodic {n x m) -matrix-valued functions M;j(x) in W^, l,j = l,...,d, 
such that 

MijeH\n), f Mi,(x)dx = 0, M,,(x) = -M,7(x), l,j = l,...,d, 

(3.2) 

d 

/Kx) = ^a,M,,(x), l = l,...,d. (3.3) 

We have 

\\Mijh,(n) ^ (2r-o)-' {Wfih.m + , l,j = 1, • • • ,d. (3.4) 

Proof. Let $i(x), / = 1, . . . ,d, be L-periodic (n x m)-matrix-valued func- 
tions in such that 

A$;(x) = /Kx), / $;(x)dx = 0, l = l,...,d. (3.5) 
Jq 

The solvabihty condition // = for the equation in (3.5) is satisfied. The 
solution exists, it is unique, and <!>/ G i/^(rj), since /; G L2(ri). We put 

M,,(x) := a,$Kx) - 9/^,(x), l,j = l,...,d. (3.6) 

Then M^^- G #^(0), 1% = 0, and M;j(x) = -Mj7(x). We have 

d d / d 

^a,M,,(x) = Y.(9]M^) - djdi^ji^)) = A$Kx) - di 5^a,$,(x) 
j=i 3=1 \i=i 

(3.7) 

Note that, by (3.1) and (3.5), 

d d / d 

= Y.d,f,i^) = 5^9,A$,(x) = A j;5,$,-(x) 

i=i i=i \i=i 

Thus, X]^=i^i^i is a periodic harmonic function with zero mean value. 
Hence, Ej=iO'j^j = 0. Now from (3.5) and (3.7) it follows that /;(x) = 
Ylj=i^j^lji^)- Relations (3.2) and (3.3) are proved. 

It remains to check estimate (3.4). Using (3.5), it is easily seen that 
satisfies the estimate 



|D$ 



i\\l2{n) ^ ll/dlL2(n)ll^dlL2(f2)- 
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Note that for any F-periodic function ip G H^{Q.) with zero mean value we 
have ||<^||L2(n) ^ (2ro)~lDy3||L2(n). Hence, 

l|D^dlL2(n) ^ (2ro)-'||/;||L2(n), l = l,...,d. 
Together with (3.6) this imphes (3.4). • 

3.2. Estimates for integrals over the neighborhood of the bound- 
ary. The following statement is rather standard (see [PSu2, Lemma 5.1]). 

Lemma 3.2. Let O dW^ he a hounded domain with the houndary of class 
C^. Denote Be = {x G C> : distjx, 90} < e}. Suppose that the numher 
£o € (0, 1] is such that B^^ can he covered hy a finite number of open sets 
admitting diffeomorphisms of class rectifying the houndary dO. Then 
for any function u G H^{0) we have 

/ \ufdyi^/3e\\u\\Hi(o)\\u\\L2{0): < e ^ eq. 

J Be 

The constant P = P{0) depends only on the domain O. 

The following statement is an analog of Lemma 2.6 from [ZhPas]; in the 
present form it was proved in [PSu2, Lemma 5.3]. 

Lemma 3.3. Let O C R'^ he a hounded domain with the houndary of class 
C^. Denote (dO)^ = {x G M'^ : dist{x, 90} < e}. Suppose that the numher 
Eq G (0, 1] is such that {dO)cQ can he covered hy a fi,nite numher of open sets 
admitting diffeomorphisms of class rectifying the houndary dO. Let S^ 
he the operator (1.17). Let /(x) he a T -periodic function in such that 
f e L2(ri). Then for any function u G ii"^(M'^; C"*) we have 

l/'^l^l'S'eUpdx ^ ,5*e|ri|"^||/|||2(f2)||u||^i(Rd.cm)||u||i2(K'*;C'")' 

< £ ^ £l, 

where £i = £o(l + ri)~^, 13^, = /3°(1 + ri), the constant depends only on 
the domain O, and, 2ri = diamO. 

In what follows, the number eq is chosen as in Lemma 3.3. Then eq also 
satisfies the assumptions of Lemma 3.2. 

3.3. The Dirichlet problem. We need the properties of solutions of the 
Dirichlct problem for the equation —As + s = f in the domain O. Recall 
that the space iJ~^(0;C'^) is defined as the space dual to ifQ(0;C") with 
respect to the L2(0; C")-duality. If f G ii"-i(0;C") and rj G i7^(0;C"), 
the symbol (f, ■'7)L2(e') is understood as the meaning of the functional f on 
the element 7). Herewith, 

l(f,^)L2(0)! ^ l|f|li/-i(Ci)ll^lbi(0)- (^•^) 
The following lemma is checked in a standard way. 

Lemma 3.4. Let f G H-^{0; C"), and let s G Hl{0; C") he the generalized 
solution of the Dirichlet prohlem 

-As + s = f inO, s\do = 0. (3.9) 
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Then the following "energy estimate" is valid: 

MIh^o) ^ l|f|l//-i(ci)- (3.10) 
Proof. By definition of the generalized solution of the problem (3.9), we 



have 

d 



J^(^^{Dis,Dirj) + {s,rj)jd^ = Jjf,rj)d^, ^rj e H^{0 (3.11) 

Substituting 77 = s in (3.11) and using (3.8), we obtain (3.10). • 

Next statement concerns the Dirichlet problem for the homogeneous equa- 
tion; it is deduced from Lemma 3.4. 

Lemma 3.5. Let 4> G H^{0;C''), and let h G H^{0;C'') be the solution of 
the problem 

-Ah + h = inO; h\do = (pldo- (3.12) 

Then we have 

\Mm{o) ^ M4>\\m{o)- (3-13) 

Proof. By (3.12), the function h - (j) e H^{0;C'') is the solution of the 
problem 

(-A + 7)(h-0) = f inO; (h- 0)1^0 = 0, 

where f = -(-A + 7)0. Since {f,v)L,{0) = -{<I>,v)hHO), V e H^{0;C^), 
then f G H-^{0; C"), and we have 

iifii \{4>-v)in(o)\ ||,|| 

¥\\h-^0) = sup ——- = 11011^1(0)- (3-14) 

By Lemma 3.4 and (3.14), 

W^ - <f>\\m{o) ^ l|f ll//-i(c>) = 
This implies (3.13). • 



§4. HOMOGENIZATION OF THE NEUMANN PROBLEM: MAIN RESULTS 

4.1. Formulations of the main results. Now we formulate main results 
of the paper. We start with approximation of the operator {AN,e + 
in the L2-operator norm for small £. 

Theorem 4.1. Suppose that O CW^ is a bounded domain with the boundary 
of class C^'^, and the matrix g{x) and DO 6(D) satisfy the assumptions of 
Subsection 1.2. Suppose that Condition 2.1 is satisfi^ed. Let A be subject to 
the restriction (2.6). Let u^ be the solution of the problem (2.8), and let Uq 
be the solution of the problem (2.28) with F G L2(C';C'*). Then there exists 
a number £1 G (0, 1] depending on the domain O and the lattice F such that 

l|U£ - U0||l2(O;C") ^ Co(A)£||F||i2(o.cn), < £ ^ £1, (4.1) 
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or, in operator terms, 
||(-4.Ar,£ + XI)"^ - {A% + XI)~'^\\l2(0;C^)^L2{0;C^) ^ Co(A)e, < e ^ £1. 

(4.2) 

The constant Co(A) depends on m, d, oq, ai, H^Hloo) \\g ^Wlocj -^j P^' 
rameters of the lattice T, the constants C\ and C2 from the inequality (2.2), 
and the domain O. 

Our second result concerns approximation of the resolvent of An,z in 
the norm of operators acting from L2(0;C") to -ff^(C';C"). In order to 
formulate this result, it is necessary to introduce the corrector. 

We fix a linear continuous extension operator 

Po : H^{0; C") ^ H'^iR'^; C"), (4.3) 
and put uq = Po^o- Then 

l|U0||//2(I[jd.cn) ^ C'c'||uo||_ff2(0.cn), (4.4) 

where Co is the norm of the operator (4.3). 

Let /Se be the smoothing Steklov operator defined by (1.17). By Rq we 
denote the operator of restriction of functions in onto the domain O. 
The corrector for the Neumann problem is given by 

KnAs) = RomSeb(D)Po{A% + (4.5) 

The operator b{'D)Pa{A% + XI)^^ is a continuous mapping of L2(C';C") 
to i?^(R°'; C™"). As was mentioned in Subsection 1.6, the operator [A'^jS'e 
is continuous from H^{W^;U^) to i7^(R'';C"). Hence, the operator (4.5) is 
continuous from L2(0;C") to iJi(C';C"). 

The first order approximation to the solution is given by 

V, = {A% + Xiy'F + eKN,x{£)F. (4.6) 

Consider the following fimction in W^: 

= So + eA' S,{b{'D)uo) = uo + eKxie)uo, (4.7) 

Here Kx{£) is the operator (1.25). Then 

= Ve\o- (4.8) 

Theorem 4.2. Suppose that the assumptions of Theorem 4.1 are satisfied. 
Let be the function defined by (4.5), (4.6). Then there exists a number 
£i G (0, 1] depending on the domain O and the kittice T such that 

W - Ve||j/i(0;C») < e{X)e^l'^\\F\\L^^^O;Cr^), < £ ^ £1, (4.9) 
or, in operator terms, 

\\{AN,e + A/)-' - {A% + AI)-1 - eKN,x{£)\\L2{0;Cr^)^m{0;Cn) ^ C{X)e^'\ 

(4.10) 

For the flux := g'^6(D)u£ we have 

lIPe - r'5e6(D)5o||L2(0;C-) ^ C'(A)£V2||F||^^(0.c„), < £ ^ £i, (4.11) 
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where ^(x) is the matrix (1.9). The constants C(X) and C'(A) depend on m, 
d, ao, ai, WqWl^, \\9~^\\Looj ^i ihe parameters of the lattice V, the constants 
C\ and C2 from the inequality (2.2), and the domain O. 

Now we distinguish the case where the correetor is equal to zero. Next 
statement follows from Theorem 4.2 and Proposition 1.2. 

Proposition 4.3. Let he the solution of the problem (2.8), and let Uq 
he the solution of the problem, (2.28). // = g, i. e., relations (1.13) are 
satisfied, then A = 0, K^^\{e) = 0, and we have 

l|U£ ~ Uo||ifi(0;C") ^ C('^)^"^^^I|F||l2(C';C")) < £ ^ £1. 

4.2. The first part of the proof: introduction of the boundary layer 
correction term. The proof of Theorems 4.1 and 4.2 relies on application 
of the results for the homogenization problem in (Theorems 1.6, 1.7, 1.8) 

and on estimation of the boundary layer correction term. 
Consider the function uq = PoUq € H'^{W^\U^). Clearly, 

F:=.4°Uo + UoGL2(M'^;C"). 

Herewith, F|o = F. Using the Fourier transformation and (1.4), (1.15), we 

obtain 

^ / ("ib^ll^l^ + l)^|uo(^)|^c?l ^ (max{ai||5(||L^,l})^ ||uo||^2(iR.dv 

(4.12) 

Here uo(^) is the Fourier-image of the function uo(x). Combining (4.12) 
and (4.4), (2.33), we arrive at 

1 1 I La (US'*) ^ ^3(A)||F||^2(0)' 

C3(A) = CxCoraa;K{ai\\g\\L^,l} . 

Let Ue G H^{W^; C") be the generalized solution of the equation 

Ae'^e + AS, = F. (4.14) 
Applying Theorems 1.6, 1.7, and 1.8, we obtain 

l|u£-uo||L2(iR'i) ^ C'i(A)£||F||z,2(]Rd), £>0, (4.15) 

- ^e\\m{W.d) ^ C'2(A)£||F||i2(iRd), £ > 0, (4.16) 

||/6(D)S, - g'Seb{T>)^4L,{^a) ^ C5(A)£||F||i^(„d), £ > 0. (4.17) 

Here is defined by (4.7). 

Consider the difference Ug — in the domain O. By (2.8) and (4.14), 
this difference is the solution of the Neumann problem for the homogeneous 
equation: 

(^, + A7)(ue-U,) = inO; 

55 (u, - S,) = -0, on dO, ^ ■ ^ 



(4.13) 
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where := df^Ue- Since Ue\o G H^{0;C'^) and AeU^ = F - Aug G 
L2{0; C"), we can apply the definition of the conormal derivative from Sub- 
section 2.4. According to (2.20) and (2.23), the element G H-^/'^{dO; C") 
is defined by 

i^s,'P)L,idO)= [ ((/6(D)u„6(D)r<^) + A(u„r¥'))dx- / {F,Tcp)d^, 
Jo Jo 

<f G H^/'^{dO-C). 

(4.19) 

Here T is the operator (2.15) defined in Subsection 2.3. 

We wish to approximate the difference — by a function Wg, which 
is solution of the Neumann problem for the homogeneous equation with a 
simpler given function in the boundary condition. The choice of is 
prompted by approximation of the solution (4.15) and approximation of the 
flux (4.17). We define p^ G H-^'^{dO;U') by the relation 

(p.,v)L,(ao)= / ((r5.6(D)So,5(D)r¥.) + A(uo,T(^)-(F,T(p))dx, 
Jo 

V? G H^/^{dO;C). 

(4.20) 

Clearly, the right-hand side of (4.20) is an antilinear continuous functional 

of G H^^^{dO;C"'), so that definition (4.20) is correct. Let us estimate 
the norm — P£\\H-^/2{dO)- From (4.19) and (4.20) it follows that 

{^,-p„^)L,i90)= [ {9'b(D)n,-rSsbiB)uo,b(D)T<p)dK 
Jo 

+ A /" (u, - uo. Tip) dx, <p G H^/\dO; C"). 
Jo 

Together with (4.15) and (4.17) this implies 
\{fl^e- Pe,'P)L2idO)\ 

^ £||Fk,(M<i) (C5(A)||6(D)r(^|U,(o) + ACi(A)||T<^|U,(o)) , (4.21) 
¥5 G iJ^/2(aC';C"), £>0. 

Combining (4.21) with (1.2), (1.5), (2.16), and (4.13), we obtain 

IKV'e - Pe,'P)L2{dO)\ ^ C4(A)e||F||i^(o)||<^||^i/2(ao), 

if e H^/^ {do £>0, 

where C^iX) = C3(A)ci (C5(A)(dai)V2 + ACi(A)). Prom here and (2.22) it 
follows that 

||Ve-P.lli/-l/2(aO) ^C4(A)£||F|U,(0), £>0. (4.22) 



28 



T. A. SUSLINA 



Now we introduce the correction term Wg G H^{0;C^) which is the gen- 
erahzed solution of the Neumann problem 

{As + AJ)w, = in O; 

= p, on dO. ^ ■ ' 



Lemma 4.4. Let Ug be the solution of the problem (2.8), and let be the 
solution of the equation (4.14). Let-w^ be the solution of the problem (4.23). 
Then for e > we have 

\\Ue - Ue + We\\m{0) ^ CbW^W^WL^iO), (4-24) 

where C5(A) = c^^ciC^^X). 

Proof. By (4.18) and (4.23), the function — + is the solution of 
the problem 

{Ae + A/)(Ue -Ue+We)=0 in O] 
df,{us - Ue + We) = p^-lj)^ on dO . 

Applying Proposition 2.4, we obtain 

||U£ - Ue+^e\\m{0) ^ (\^Cl\\Pe " V'£llif-i/2(90)- 

Together with (4.22) this implies (4.24). • 

Conclusions. 1. Combining (4.8), (4.13), (4.16), and (4.24), we arrive 

at 

||u,-v, + w,||^i(e,) ^ (C5(A) + C2(A)C3(A))£||F||i^(0), £ > 0. (4.25) 

Hence, it is clear that, in order to prove main estimate (4.9) from Theorem 
4.2, it suffices to estimate \\^E\\m{0) C^^^'^\\'^\\l2{0)- 

2. Replacing the H'^{0\ C")-norm by the L2(P; C")-norm in the left-hand 
side of (4.24), we obtain: 

||Ue - + We||L2(0) ^ C5(A)£||F||i2(o), £ > 0. (4.26) 

Combining (4.26) and (4.13), (4.15), we arrive at 

||u, - uo + w,|U,(o) ^ (C5(A) + Ci(A)C3(A)) £||F|U,(o), e > 0. (4.27) 

This shows that, in order to prove estimate (4.1) from Theorem 4.1, it 
suffices to estimate ||we 11^2(0) terms of C£||F||^2(e')- 

Thus, the problem is reduced to estimation of the function which can 
be interpreted as the boundary layer correction term. First wc estimate the 
norm of in H^{0;C^), and so prove Theorem 4.2 (sec §5). Next, using 
already proved Theorem 4.2, we will estimate the norm of in L2(0;C") 
(see §6). 
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§5. Proof of Theorem 4.2 

5.1. According to Subsection 2.4, the generalized solution G i/^(0;C") 
of the problem (4.23) satisfies the identity 

J ((/fe(D)w„6(D)r?) + A(w„77))dx = (p„7T?)i,(ao), & H\0;Cn. 

(5.1) 

Using (4.20) and (2.30), we represent the right-hand side of (5.1) as 

{pe,lv)Mao) = [ {{rSeb{D)uo,b{D)T^ri) + X{uo,Tjri) - {F,Tjri)) dx 
Jo 

= [ {rS,b(D)uo-g%(D)no,b(D)T^rj). 
Jo 

(5.2) 

Prom (5.1) and (5.2) it follows that 

y"((/fe(D)w„6(D)r?) + A(w„r?)) dx = I,[r?], ^rj G H\0;C'), (5.3) 
o 

where we denote 

Uv]-= I (r^sKD)So-A(D)uo,6(D)r777)dx, 77 G //^(O; C"). (5.4) 
Jo 

For further estimates, it is convenient to represent the functional (5.4) as 

X,M=X«M+Xf)[.7], (5.5) 

where 

^i'^ H = I (/5.fe(D)uo - 5°KD)uo, h{T>)T-fini) dx, (5.6) 
Jo 

4^\v\ = I {{r-g°)Seb{-D)^,b{D)Tjrj)d^. (5.7) 
Jo 

The term (5.6) can be easily estimated with the help of Proposition 1.4 
and relations (1.2), (1.5), and (1.15): 

\4'^[V]\ ^ b°|||(5.-/)6(D)uo||i,(M.)||6(D)r7r7||i,(o) 

^ 49\\L^n\\m(D)uo\\L,^^,)idaif/^\\T^T]\\Hi^O), £ > 0. 
By (1.4), (2.33), and (4.4), 

^ 1/2 

||D6(D)uo||L2(Md) l|uo||/f2(Rd) 

^ al^^Co ||uo||i/2(o) ^ ai^^C'o^A||F||L2(0). 
Now, combining (5.8), (5.9), (2.13), and (2.16), we arrive at 

\iP[v]\^CeiX)8\\F\\L,^0)Mm(0), V&H\0;C-), e > 0, 
Cq{X) = d^/^riai\\g\\L^CoCxCiCi. 



(5.8) 



(5.9) 



(5.10) 
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5.2. Analysis of the term . Transform the term (5.7), using (1.2): 

„ d 

^i'^W= / ^(/f.5,6(D)5o,Ar7r7)rfx, (5.11) 

where we denote 

/Kx):=6r(5(x)-/) = 6;(5(x)(6(D)A(x) + l^)-5°), Z = l,...,d. (5.12) 

Then //(x) are F-periodic (n x m)-matrix- valued functions, fi G i^2(^^); and, 
by the definition (1.8) of the effective matrix, we have /; = 0. Finally, from 
the equation (1.7) for A it follows that 

d 

Y,Difi{^) = 6(D)*3(x)(6(D)A(x) + 1„) = 0. 
1=1 

Thus, the set of functions //(x), / = l,...,(i, satisfy the assumptions of 
Lemma 3.1. By this lemma, there exist matrix- valued functions M/j(x), 
Z, j = 1, . . . , d, such that relations (3.2), (3.3) are valid. Then 

d 

/f(x) = £5^a,Mf,.(x). (5.13) 

We will need estimates for the L2(rj)-norms of the functions (5.12). From 
(1.5), taking into account that ^ — ^f'^ is the orthogonal projection of ^ onto 
the space orthogonal to constants, we obtain 

II ^ II / 1/2||~ On ^ 1/2||~|| 

\\ji\\L2m^a^ Il5 - 5 llL2(n) ^ «! llffllL2(n)- 
Combining this with (1.2), (1.5), (1.9), and (1.10), we have 

\\fi\\L,(,n) ^ cc]'''\\9\\lM^''' + (ci«i)'/'l|DA|U,(f^)) ^ IJ^I^'C, Z = 1, . . . , d, 

(5.14) 

where £ = Q:J^^||5f||L^ (1 + ((im)"^/^Q;}^^aQ ||5f~"^||^'^^). Hence, by 

(3.4), 

^ro-^|J^|V2c, Z,j = l,...,d. (5.15) 

By (5.13), 

d 

/f5,6(D)uo = (a,Mf.) S,6(D)uo 

(5.16) 

d d ^ ' 

= eY,dj (M,^.5,6(D)uo) - e ^ Mf^. 5, 6(D) 5,- uq. 
j=i j=i 

According to (5.16), the term (5.11) can be represented as 

lP[rj]=iP[rj]+lP[rj], (5.17) 
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where 



^i'^ [V]=e [ Yl i^j {M[jSeb{T>)no) , DiTjrj) dx, (5.18) 

„ d 

W = / E {Mf^S,h{Ti)djU^, DiTjrj) dx. (5.19) 

The term (5.19) can be easily estimated by using Proposition 1.5: 
d 

I^PNI ^ e ^ |J^rV2||M,,-|U,(f,)||6(D)5,-uo||i,(K.)||Ar7»7llL,(0)- 

Together with (2.13), (2.16), (5.9), and (5.15) this imphes: 

\WM^C7{X)s\\F\\L^^O)MmiO), V&H\0;C-), e > 0, 

I/O ( o.zU I 

C7(A) = dro'<tay CoCxCici. 

~(2) 

5.3. Estimate for the term . Recall the notation {dO)^ for the e- 
neighborhood of the boundary dO in W^. Assume that < e ^ £i, where 
the number ei was defined in Lemma 3.3. We fix a smooth cut-off function 
Oeix) in M'^ such that 

0,eC^{R''), suppOeCidO),, 0^^,(x)^l, ^^^^^ 

6'e(x)|ac' = 1, e|V6'£(x)| ^ k = const. 

Obviously, such a function exists. The constant k depends only on the 
domain O. 
Note that 

d 

^ {dj (M;^.(l - e,)S,b(D)uo) , DiTjrj) dx = 0, r? G H\0; C"). 



(5.22) 

Indeed, since the left-hand side of (5.22) is an antilinear continuous func- 
tional of r? G i?^(0;C"), it suffices to prove this identity for some dense 
set in H^{0;C'^). Assume that r] G H'^{0;C^), and, integrating by parts, 
rewrite the left-hand side of (5.22) as 

„ d 

-i / ((1 - e,)Seb{B)uo, V {M[^rdjdiT^ri) dx. 

This expression is equal to zero since Ylf j=ii^ij)* ^j^^l = ^ due to Mij = 
-Mji (see (3.2)). 

By (5.22), the term (5.18) takes the form 

„ d 

H=e E (^.Mf^-5sKD)So) , DiT^r^) dx. (5.23) 

•'''1,3=1 
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Consider the following expression 

d d 

d d 

+ e Y,(.dje,)MfjS,b{^)uo + Y,idjMijYS,b{^)uo. 
Take into account that Y!j=idjMij = fi (see (3.3)). Then 

^j('\s) + jP{s) + 4'\s), (5.24) 



■13= 

d 



L2{0) 



where 

d 

4'\s) =eJ2P,mjSM'D)djUo)h,iO), (5.25) 

d 

J^is) = ||(a,e,)Mf^.5,6(D)So||L,(o), (5.26) 
j=i 

4^\e) = \\ej!S,b{^)uoh,io)- (5.27) 
In order to estimate the term (5.25), apply (5.21), Proposition 1.5, and 
relations (5.9), (5.15): 
d 

J«(£) ^ e^|J^|-V2||M,,-||^,(o)||fe(D)5,-uo|L,(M.) ^ £«(A)£||F||^,(o), 

(5.28) 

where e:(i)(A) = ro^dd^/^al^^CoCx- 

The term (5.26) is estimated with the help of (5.21) and Lemma 3.3: 

{4'\e)f^d^y' [ \Mi/\S,b{U)uo\^d^ 

(5.29) 

^ £ci«:2/3,|n|-i^ ||M,,-||i^(^)||S,6(D)uo||l,i(Kd), < £ ^ £1. 

From (1.19), (1.4), (2.33), and (4.4) it follows that 

\\Seb{Y))uo\\Hi{ni) ^ l|fc(D)uolbi(Md) ^^^^^ 

^ «i'^^l|So|lH2(]Rd) ^ '^^i'^CoCx\\F\\l2{0)- 
Combining (5.29), (5.30), and (5.15), we obtain 

j(')(£)^£(2)(^)^l/2||p||^^^^^^ 0<£^£l, (5.31) 

where = KdpV'^rQ^(La^'^CoCx. 



HOMOGENIZATION OF THE NEUMANN PROBLEM FOR ELLIPTIC SYSTEMS 33 

The term (5.27) is estimated in a similar way by using (5.21) and Lemma 
3.3: 

(jf (e))2^ / |/f|2|5,6(D)So|^c/x 

J{dO), (5.32) 

^ £/3.|J^|-i/i||i^(^^||5,6(D)uo|||i(ud), < £ ^ £1. 
Prom (5.14), (5.30), and (5.32) it follows that 

4^\e) ^ £(^)(A)£V2||F||^^(o), < £ ^ £i, (5.33) 

where (t^^\X) = pl^^fta/^CoCx- 

Finally, relations (5.24), (5.28), (5.31), and (5.33) imply that 

d , _ 

^ £(A)£V2||F||i2(0), 0<£^£i, 

L2{0) 

(5.34) 

where C(A) = £(^(A) + (L^'^\X) + (t^^\X). Combining (2.13), (2.16), (5.23), 
and (5.34), we arrive at 

\iP[v]\ ^ C8(A)£V2||F||i^(o)||r;||^i(o), G H\0;C''), < £ ^ £i, 

(5.35) 

where CsiX) = (t{X)d^/^ciCi. 

5.4. Completion of the proof of Theorem 4.2. Now relations (5.5), 
(5.10), (5.17), (5.20), and (5.35) imply that 

\Urj]\^Cg{X)e'/''\\FU,^0)\\v\\miO), V ^ H\0;C^), < £ ^ £i, 

(5.36) 

where C9(A) = C6(A) + CjiX) + Cs{X). 

From (5.36) and (5.3) with rj = it follows that 

J ((/6(D)w„6(D)we) + A|we|2) dx^C9(A)£V2||F||i^(o)||w,||j^i(o) 
o 

for < £ ^ £1. Together with (2.7) this yields the desired estimate 

Mm(o) < c-^CgiXy/^\\F\\L,^a), < £ ^ £i. (5.37) 

As a result, (4.25) and (5.37) imply (4.9) with the constant 

C(A) = C5(A) + C2(A)C3(A) + c^'Cg{X). 

It remains to obtain approximation (4.11) for the flux. By (4.9), taking 
(1.2) and (1.5) into account, wc have 

Hp, -/6(D)v,||i,(0) ^ \\g\\L^{daif/^C{X)e^/^\\F\\L,^o), < £ ^ £i. 

(5.38) 

Let V£ be defined by (4.7). Note that, by (4.8), 
||/6(D)v, - rS.6(D)5o|U,(o) ^ ||/6(D)v, - S,b{-D)uo\\L,(Rdy (5.39) 
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Let US use the proof of Theorem 1.8. By (1.47), 

||/6(D)v, - 5^5,6(D)5o||i,(Kd) ^ C6e||D6(D)5o||i,(K<i), £ > 0. (5.40) 
Prom (5.39), (5.40), and (5.9) we obtain 

||/6(D)v, - g'S,b{-D)uo\\L,(0) ^ Cea\/^CoCxe\\F\\L,^o), e > 0. (5.41) 

Finally, (5.38) and (5.41) imply the required estimate (4.11) with the con- 
stant C'(A) = \\g\\L^{dai)^/'^C{X) + CeaY^CoCx. This completes the proof 



of Theorem 4.2. • 



56. Proof of Theorem 4.1 



6.1. As was mentioned in the end of §4, for the proof of Theorem 4.1 we 
need to estimate the norm ||w£||j^2(0) by C£||F||£,2(C')- Comparing (5.5), 
(5.10), (5.17), and (5.20), we see that 

|I,[r7]| ^ (C6(A) +C7(A))£||F|U,(o)||r7||^i(o) + \iP[ri]\, 

rjeH^iO-C), £>0. 

In the identity (5.3) which is valid for any function rj G //^(C';C"), now 
we choose 77 in a special way. Namely, let $ G L2(C';C"), and let r}^ G 
H^{0\C^) be the generalized solution of the Neumann problem 

5(D)*/KD)T7e + At7, = * in O; 0^^7^190 = 0, (6.2) 

i. e., T), = {AN,e + A/)-i*. By (2.11), 

\\'ne\\m{0)^cl^\\^\\L^(0)- (6-3) 

By definition of the generalized solution of the problem (6.2), we have 

/ ((/6(D)w„6(D)r7,) + A(w„r7,)) dx= /(w„*)dx. (6.4) 
Jo Jo 

Combining (6.4) and (5.3), we obtain 

(w„*)dx = X,[r/,], *GL2(C»;C"). (6.5) 



10 

Now (6.1), (6.3), and (6.5) imply that 

|(w.,*)l,(0)| ^ {C^{X)+C,{X))clh\\F\\L,(o)m\L,(0) + \W\Ve\\, 

* G L2(0;C"), £ > 0. 

6.2. Analysis of the term X^\r}i,]. Clearly, the required estimate for 
1 1 We 1 1 £^2(0) '^ill be proved, if we succeed to estimate the last term in the 
right-hand side of (6.6) by Ce\\F\\L^^o)\\^\\L^^oy 

We apply the (already proved) Theorem 4.2 in order to approximate the 
solution -q^ of the problem (6.2) in the iJ^(C); C")-norm. Denote by t^q the 
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solution of the corresponding homogenized problem: ryg = (.4^ + A/)^^$. 
Let T7o = PoVo- Note that, by (2.34), 

ll^olli/2(Kd) ^ Co 1 1 r7o I Iff 2(0) ^ CoCx\\^\\l2{o)- (6-7) 
Prom (4.10) it follows that 

||77, -rjo- eA'S,b(D)?iJm^o) ^ CiX)e^/^\ML,(o), < e ^ si. (6.8) 
Represent [^7^] (see (5.23)) as the sum of three terms 

iPiVe] = iP[{AN,e + A/)-^*] = + + jPm, (6.9) 

where 

„ d 

■^^ i,j=i 

(6.10) 

„ d 

JP [*] = ^ / E {0,M[^Seb{D)uo) , Ar777o) d^, (6.11) 
Jo 



d 

lO 



„ d 

jPI^] = £ / E {0eM[jSeb{Y>)no),DiT-fieA'S,b{n)7)o)) dx. (6.12) 



The term (6.10) is easily estimated by using (6.8), (2.13), (2.16), and 
(5.34): 

<cW(A)e||F||i,(0)||*||^,(0), *GL2(0;C"), < e ^ si, 

(6.13) 

where C(i)(A) = e:(A)di/2?iciC(A). 

In order to estimate the term (6.11), we apply (5.21) and (5.34): 

d 

^ e^/^5(A)||F||i2(C') 

1=1 

According to Lemma 3.2, 

W^lTlVoWl^iB,) ^ P4^lT'yVo\\mio)\\DiT'yr]o\\L2{0), < e ^ £o- 
Hence, 

d 

E WDiTiVoWmb,) ^ £'/'(^^i)'/'r7^ollJ^2(o)||T7r?oll^/i(o), < e ^ sq. 
1=1 
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Applying estimates (2.13), (2.16), (2.18), (2.32), and (2.34), we obtain 
d 

^V2||„ ||V2 ||„ ||l/2 

1=1 



Yl llAr7T7ollL2(i3e) ^ £^^^(/3t^)^^^(c2C2CiCi)^/^||77o||^2(o)||r7olli^i%) 

(6.15) 



Now (6.14) and (6.15) imply that 

\jPm^C^^\X)e\\F\\L^^o)\mL,io), *GL2(0;C"), < £ ^ eo, 

(6.16) 

6.3. Estimate of the term It remains to estimate the term 

(6.12). Denote 

z, := TjieA'S,b(D)?io). (6.17) 
By (5.34), the term (6.12) satisfies 

|ji3)[*]K£V2£(A)||F||^^(0)^||Az.|U,(0), 0<£^£i. (6.18) 

1=1 

By definitions of the operators 7 and T (see Subsection 2.3), the function 
(6.17) is the generalized solution of the Dirichlet problem 

- Azs + Ze = in O; 

ZeldO = eA^Seb(D)7io\9o. ^^'"^^^ 

Now we estimate the norm of z^ in i^^(0;C"'). Consider the following 
function in W^: 

0,(x) = £0,(x)A^(x)(Se6(D)?7o)(x). (6.20) 

Taking (5.21) into account, we can rewrite the problem (6.19) in a different 
way: — Az^ + z^ = in O and z^Iqq = (p^ldo- By Lemma 3.5, 

heWnHo) ^'^HeWmHo)- (6-21) 
Thus, the question is reduced to the following statement. 
Lemma 6.1. Let 0^ be the function (6.20). Then we have 

\\<t>e\\m{0)^Cio{X)e'/^\ML,^0), 0<£^£i. (6.22) 

Proof. The proof is similar to that of Lemma 7.4 from [PSu2]. 

First, we estimate the L2(C'; C")-norm of the function (6.20). By (5.21) 
and Proposition 1.5, we have 

Ue\\L2{0) ^ £||A"5',5(D)^7o||i2(K<i) ^ £\^\~^^^\\M\L2{n)\\H^)Vo\\L2{R'iy 

(6.23) 

Prom (1.4) and (6.7) it follows that 

IIKD)T7ollL2(]Rd) ^ a}^^ll^ollifi(K'') ^ (h.^^CoCx\\^\\L2{0)- (6-24) 
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Now relations (6.23), (6.24), and (1.11) imply that 

12(0) «^ 7o(A)e^||*||i2(C,), (6.25) 

where 70(A) = m{2ror^aia^%\\Lj\9-'\\L^C^oCl 
Consider the derivatives 

+ e0eA'Seb(D)dj7io, j = l,...,d. 

Then 

l|D0sllL(o) ^ 3e2^|Ve,|2|A^5,6(D)?7o|2dx 

+ 3/ \{T>Ay\^\9,Seb{U)?]o\^d^ + 3e^y2 [ \e,\^\A'S,b{n)Dj7jQ\^ d^. 
Jo j^i-^^ 

(6.26) 

Denote the summands in the right-hand side of (6.26) by 2li(e), 2t2(£), and 
213(e), respectively. 

The easiest task is to estimate 213(e). By (5.21) and Proposition 1.5, 

d 

213(e) ^ 3e2^ ||A^5,6(D)D,?7ollL(Mc*) 
i=i 

^3e2|^7|-i||A||i^(^)||D6(D)?7of 
Similarly to (5.9), by (6.7), 



(6.27) 



1^2(0)- 

Prom (6.27), (6.28), and (1.11) it follows that 

2l3(£)^73(A)e2||*||i^(0), (6.29) 

where 73(A) = 3m{2ro)-^a^a^^g\\Lj\9-'\\L^ClCl 

The first term in the right-hand side of (6.26) can be estimated by using 
(5.21) and Lemma 3.3. For < e ^ ei we have 

2li(e) ^ 3k2 / \A' S,b(D)Jjof dx 

J{dO), (6.30) 

^3/c2/3,£|J7|-iA||2^(^)||6(D)^7ol|^i(M.)- 

Note that, by (1.4) and (6.7), 

\\K^)riQ\\m(9.i) ^ a}^^ll^7oll//2(Rd) ^ (A'^CoCx\\^\\l2{0)- (6-31) 
As a result, (6.30), (6.31), and (1.11) yield 

2ti(e) ^7i(A)e||*||i^(0), < e ^ ei, (6.32) 

where 71(A) = ^k-" ^m{2r^)-^aia-^^\\g\\Lj\9-^\\L^ClCl. 
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It remains to consider the second term in the right-hand side of (6.26). 
By (5.21), 

212(e) ^ 3 [ \{T>Arf\S,biB)7jofd^. 

J{dO)e 

Applying Lemma 3.3, for any < £ ^ £1 we obtain 

2l2(£) ^ 3^.£|0|-i||DA||i^(^)||6(D)r7o||^i(K<i). 
Together with (1.10) and (6.31) this impUes that 

2t2(e) ^72(A)£||*||i^(o), 0<£^£i, (6.33) 

where 72(A) = 3P*rnaia^^\\g\\L^\\g-'^\\L^Cl,Cl. 

Finally, from (6.26), (6.29), (6.32), and (6.33) it follows that 

\\^<t>e\\l,(o) < 2ti(e) + 2l2(£) + 2t3(e) 

^ (71(A) + 72(A) +73(A))£||*||i^(0), < £ ^ £1. 
Together with (6.25) this yields (6.22) with the constant 

Cio(A) = (70(A) + 71(A) + 72(A) + 73(A))'/' . . 

From (6.21) and Lemma 6.1 it follows that 
^2Cio(A)£V2 11*11 

^2(0)' < £ ^ £1. (6.34) 

Now (6.18) and (6.34) imply that 

\jPm^C^^\X)s\\FU^^o)\mL,ia), *GL2(0;C"), < £ ^ £1, 

(6.35) 

where C(^\X) = 2€{X)d^/^Cio{X). 

6.4. Completion of the proof of Theorem 4.1. Let us summarize the 
results. From (6.9), (6.13), (6.16), and (6.35) it follows that 

\iP[Ve]\ ^ (CW(A) +C(2)(A) +C(=')(A))£||F|U,(o)||*||i,(o), 

* G L2(C>;C''), 0<£^£i. 

Together with (6.6) this yields 

l(We,*)L2(0)KCii(A)£||F||^,(0)||*|U,(0), *GL2(0;C"), 0<£^£i, 

where Cii(A) = {Cg{X) +C7{X))c^^ +C^^\X) +C^'^\X) +C^^\X). The desired 
estimate follows: 



L,iO)^Cn{X)e\\F\\L,io), < £ < £1. (6.36) 



Finally, relations (4.27) and (6.36) imply (4.1) with the constant Co(A) 
C5(A) + Ci(A)C3(A) + Cii(A). Theorem 4.1 is proved. • 
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§7. HOMOGENIZATION OF THE NEUMANN PROBLEM IN THE CASE A G Lqo 

7.1. As for the homogenization problem in M"* (see Subsection 1.7), the 
smoothing operator in the corrector (4.5) can be removed under addi- 
tional assumptions on the matrix A(x). 

Assume that Condition 1.10 is satisfied, i. e., A G L^o- Denote 

K%^,{e) = mb{B){A% + XI)-\ (7.1) 

By (2.34), the operator 6(D) (^^ + A/)~^ is a continuous mapping of 
L2(0;C'") into H^{0;C"'). Under Condition 1.10 the operator [A=] is con- 
tinuous from H^{0;C"^) to H^{0;C"-) which easily follows from Propo- 
sition 1.11 and existence of a linear continuous extension operator from 
Fi(0;C™) to Fi(M"';C'"). Hence, the corrector (7.1) is a continuous map- 
ping of L2(0;C") to H\0;C''). 

Instead of (4.6), we consider another approximation of the solution of 
the problem (2.8): 

= {A% + A7)-^F + sK%^^{s)F = uo + eA^6(D)uo. (7.2) 

Theorem 7.1. Suppose that the assumptions of Theorem 4.1 and Condition 
1.10 are satisfied. Let be defined by (7.2). Then there exists a number 
£\ G (0, 1] depending on the domain O and the lattice T such that 

l|Ue - Vgllj^i^Q.Cn) ^ C(A)£^/^||F||^2(Cl;C")) 0<£^£i, (7.3) 

or, in operator terms, 

\\{AN,e + XI)-' - {A% + XI)-' - £i^^,A(£)llL,(0;C")-.i/i(0;C") ^ ^(A)^^^ 

For the flux = g^b(D)ue we have 

Hp, - r&(D)U0||L2(O;C-) ^ C'(A)£V2||F||^^((p.c„), < £ ^ £i, (7.4) 

where g'(x) is the matrix (1.9). The constants C{X) and C'(A) depend on m, 
d, ao, ai, IbllLoo) llfi'~^llioo' A, the parameters of the lattice F, the constants 
Ci and C2 from the inequality (2.2), the norm ||A||l^, and the domain O. 

Proof. According to (4.5), (4.6), and (7.2), in the domain O we have 

ve-v, = £A=(I-5,)6(D)Go. (7.5) 

Let us estimate the norm of the function (7.5) in if^(0;C"). Since the 
right-hand side of (7.5) is defined in the whole W^, it suffices to estimate its 
norm in ii"i(M'^; C"). We apply the proof of Theorem 1.12. By (1.52) and 
(1.58), 

£||A^(/-5e)6(D)uo||i,(Md) ^ 2£||A||i^||6(D)Go||L2(Md), (7-6) 

d 

^||aK£A-(/-5,)6(D)5o)||i^(«.) 

1=1 \'-') 
^ 2e' (4(1 + /32)||A||i^ +P^rl) ||D6(D)5o||i,(^.). 



1=1 
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Together with (5.9), (5.30), and (7.5) this impHes that 

\\ve - ^cWh^o) = ell^^U - Se)bi'D)uo\\Hi(Rd-) ^ Ci2(A)e||F||^2(o)> e > 0, 

(7.8) 

where C^iX) = a^^Codx (2(3 + 2/32)^/2 ||A||l^ + (2/3i)i/Vi). 

As a result, combining (4.9) and (7.8), we arrive at the estimate (7.3) 
with the constant C(A) = C(A) + Ci2(A). 

It remains to check (7.4). Prom (7.3), taking (1.2) and (1.5) into account, 
we obtain 

Hp, -/6(D)v,|U,(0) ^ e'/^dy^al/^\\g\\L^CiX)\\F\\L,^o), < e ^ e^. 

(7.9) 

By (1.2), (1.9), and (7.2), 

d 

/6(D)v, = r&(D)uo + e ^ /6,A^6(D) Auo- (7.10) 

1=1 

Using Condition 1.10 and relations (1.2), (1.5), and (2.33), we estimate the 
norm of the second term: 

d 

^ eailbllLooll^lUoo \\DjDiUo\\L2{0) 
< sdai\\g\\L^\\A\\L^Cx\\F\\L^^oy 

(7.11) 

Now (7.9)-(7.11) imply estimate (7.4) with the constant C'{X) = 

d^/'ay^\\gUj{X) + daMU\A\\L^Cx. . 

Recall that some sufficient conditions which ensure Condition 1.10 are 
given in Proposition 1.13. 

7.2. Special case. We distinguish the special case where 5° = g. Then 
5r(x) = g^ = g and, by Proposition 1.13, Condition 1.10 is satisfied. Applying 
the statement of Theorem 7.1 concerning the flux, we arrive at the following 

statement. 

Proposition 7.2. Let be the solution of the problem (2.8), and let uq be 
the solution of the problem (2.28). Let Pe = g^b(D)ue. Assume that g^ = g, 
i. e., relations (1.14) are satisfied. Then we have 

lIPe- A(D)U0||l2(O;C-) ^C'(A)£'/'||F||i,(O.Cn), 

§8. Approximation of solutions in a strictly interior subdomain 

8.1. Let O' be a strictly interior subdomain of the domain O. Applying 
Theorem 4.1 and the results for the problem in (Theorems 1.7 and 1.8), 

it is easy to obtain a sharp order error estimate for approximation of the 
solution of the problem (2.8) in if^(C'';C"). 

Theorem 8.1. Suppose that the assumptions of Theorem 4.2 are sat- 
isfied. Let O' be a strictly interior subdomain of the domain O. Let 
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5 := dist {O' ; dO} . Then there exists a number e\ £ (0,1] depending on 
the domain O and the lattice F such that 

^ £5(A)£:||F||i2(e).cn), < £ ^ £i, (8.1) 

or, in operator terms, 

\\{AN,e + A/)~^ - {A% + XI)-^ - £i^iV,A(£)||L2(0;C")^i/i(0';C") ^ <^5{X)e. 

(8.2) 

For the flux Pe = (/^6(D)ue we have 

\\Ps - 5''5e6(D)5o||i2(ci/.c™) ^ (2:5(A)e||F||i2(0;C")> < e ^ £i, (8.3) 

where g(x) is the matrix (1.9). The constants in estimates are given by 
(tsiX) = Ci(A)ri + C2(A), (t'giX) = e:'i(A)ri + (t'^iX), where Ci(A), e:2(A), 
C'i(A), and^2{X) dependonm, d, ao,ai, H^Hloo) A, i/ie parameters 

of the lattice T, the constants C\ and C2 from (2.2), and the domain O. 

Proof. We fix a smooth cut-off function C(x) such that 

CeC^iO), o^C(x)^i, 

C(x) = 1 for X G C; |VC(x)| ^ k'J'^ 



Here the constant n' depends only on the domain O. Let Ug be the solution 
of the problem (2.8), and let Ug be the solution of the equation (4.14). Then 
{Ae + XI) — Ue) = in the domain O. Hence, 



((/6(D)(u, - 5,), 6(D)77) + A(u, - S„ rj)) dx = 0, Vt? G H^{0; C"). 

(8.5) 

Substitute rj = (^(u^ - u^) in (8.5). By (1.2), 

d 

6(D) (C'(u, - u,)) = C6(D) (C(u, - u,)) + J2 biiDiOCiu, - u,), (8.6) 

1=1 

d 

CKD)(u. - S,) = 6(D) (C(u, - S,)) - ^ 6KAC)(u. - u,). (8.7) 

1=1 

Prom (8.5)-(8.7) it follows that 

J(£) := / ((/6(D) (C(u, - 5,)) , 6(D) (C(u, - S,))) + ACV. - 2,^) dx 
Jo 

= Ji{e) + J2{e), 

(8.8) 



where 



P d 

Ji(£) = - / (/6(D)(C(u,-2,)),V6KAC)(ue-3,))dx 

„ d d 

+ / {g'Y.^j{DjC,){ne-n,),Y,bi{DiC){n,-n,))d^, 



3=1 1=1 
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» d 

M^) = ia'T, &KAC)(u. - u,), 6(D) (C(u, - u,))) dx. (8.10) 
The term (8.10) is estimated with the help of (1.5) and (8.4): 

\J2{e)\ ^ ||(/)^/'6(D) (C(u, - S,)) |U,(o)||(/)^/2 J2 HDiOi^e - G.)|U.(0) 

1=1 

^ + \\9\\L^aidiK')H-''\\u, - U,||i^(o). 

(8.11) 

Similarly for the term (8.9) we obtain: 

\J,{s)\ ^ ^J{e) + 2\\gU^aid{K'f5-'\\n,-u4l^^oy (8.12) 
Now (8.8), (8.11), and (8.12) imply that 

f-2||„ ~ I, 

lL2(0)- 



J{e) ^ 6||5||L^aid(K')'<5-2||u, - Ci.||Lfnv (8.13) 



Extending the function (f^ := (^(u^ — u^) by zero to R'^ \ O, note that 
J(e) = a^lip^jcp^] + A||<^£||^^|-jgd) and apply the lower inequality (1.6). We 
obtain 

J{e) ^ ao||5"izil|D¥'£llL(Md) + ^ll¥'£llL(Md)' 

whence 

||C(u,-Se)||^i((P) ^max{aoi5-'lUoo,A-^}J(£). (8.14) 
By (4.26) and (6.36), 

||u, -u,||i,(0) ^ (C5(A)+Cii(A))£||F||i,(o), 0<£^£i. (8.15) 
Combining (8.13)-(8.15) and (8.4), we arrive at 

\\Ue - UeWn^^O') ^ IIC(U£ - ^e)\\H^O) ^ ^l(A)(^~^£||F||i2(0), < £ ^ Ei, 

(8.16) 

where 

e:i(A)=max{ao^/'|b-i||i/^,A-V2}||^||V^(6„^^)i/2^/(C5(A)+Cn(A)). 

(8.17) 

Taking (4.8), (4.13), and (4.16) into account, we check that 

115. - Vellifi(o) ^ C2(A)C3(A)e||F|U,(o), e > 0. (8.18) 

Now relations (8.16) and (8.18) imply (8.1) with the constant ^^(A) = 
e:i(A)5-i + (t2(A), where e:i(A) is defined by (8.17), and e:2(A) = C2(A)C3(A). 
It remains to check (8.3). From (8.16) and (1.2), (1.5) it follows that 

Hp, - /6(D)G,|U,(o,) ^ \\g\\L^afd'/^<ti{\)6-'e\\F\\L,^o), < e ^ Si. 

(8.19) 

By (4.13) and (4.17), we obtain 

||/6(D)S, - rSeb(D)uo\\L,io) ^ C5(A)C3(A)£||F|U,(o), 8 > 0. (8.20) 
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Now relations (8.19) and (8.20) imply (8.3) with the constant (tg{\) = 
£i(A)(5-i + £'2(A), where 

£'i(A) = |bl|L^a}/V/2(j:i(A), £^(A) = C5(A)C3(A). . (8.21) 

8.2. The case where A G L^. Similarly, under Condition 1.10, Theorems 
4.1 and 1.12 imply the following result. 

Theorem 8.2. Suppose that the assumptions of Theorem 7.1 are satisfied. 
Let O' be a strictly interior subdomain of the dom,ain O. Denote 6 := 
dist {C; (90}. Then there exists a number £i G (0,1] depending on the 
domain O and the lattice T such that 

- V£||/fi(C)';C") ^ ^5(A)£||F||i2(o.cn), < £ ^ £i, (8.22) 
or, in operator terms, 

\\{AN,e + A/)-' - {A% + A7)-l - £if^,A(£)llL2(0;C")->i/i(0';C") ^ ^sW^- 

For the flux = g^b(D)ue we have 

lIPe -3^6(D)uo||l2(0';C'") ^ ^'sW4'P\\L2iO;C"), < £ ^ £i, (8.23) 

where ^(x) is the matrix (1.9). The constants in estimates are given by 
C5(A) = e:i(A)<5-i + £2(A), £'5(A) = e:'i(A),5-i + £^(A), where £i(A), £2(A), 
(t'i{X), and€^2{X) depend on rn, d, ao, ai, HsHloo) II5~"'^IU=o> A, the parameters 
of the lattice F, the constants Ci and C2 from (2.2), the norm ||A||x,^, and 
the domain O. 

Proof. The inequality (8.16) remains true. Instead of (8.18), now we use 
the next estimate which follows from Theorem 1.12 and (4.13): 

- VelUi(o) ^ C7{\)e\\F\\L2{Rd) ^ C7(A)C3(A)£||F||i^(o), e > 0. (8.24) 

Then (8.22) is a direct consequence of (8.16) and (8.24), herewith, ^^(A) = 
(!:i(A)5-i + £2(A), where £i(A) is defined by (8.17) and £2(A) = C7(A)C3(A). 

In order to prove (8.23), we use (8.19) and the next estimate which follows 
from Theorem 1.12 and (4.13): 

||/6(D)S, -rKD)uo||L,(o) ^ C8(A)e||F||^^(i,.) 

^Cs{XmX)e\\F\\L,io), e>0. 

As a result, we obtain (8.23) with the constant £^(A) = e:'i(A)^~"^ + ^2(A), 
where e:;(A) is defined by (8.21) and £^(A) = C8(A)C3(A). • 

Remark 8.3. Theorems 8.1 and 8.2 can be applied also in the case where S 
depends on e and a strictly interior subdomain O'^ approaches the boundary 
do. For instance, if 5(e) = 0(£"), then the error estimate in (8.1) or in 
(8.22) is of order 0(£^~"). Clearly, in the case where a < 1/2 the error 
estimate in a strictly interior subdomain O'^ is better in order than in the 
whole domain O. If a ^ 1/2, then Theorems 4.2 and 7.1 give better results 
(in the whole domain O) and it does not make sense to apply Theorems 8.1 
and 8.2. 
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§9. HOMOGENIZATION OF THE NEUMANN PROBLEM IN THE CASE A = 

9.1. Statement of the problem. For applications, it is interesting to 
study the Neumann problem for the equation AN,e'^£ = F, i. e., in the 
case where A = 0. In this case, the right-hand side of the equation must 
be subject to the solvability conditions, and the solution is subject to the 
additional condition of orthogonality to the kernel of the operator. In the 
present section, we deduce the results for this case from the results of §4 
and §7. 
We denote 

Z = {ze H\0; C") : b(D)z = 0}. (9.1) 

Then Z is a (closed) subspace of H^{0;C^). Note that for sure Z contains 
the n-dimensional subspace {u € H^{0; C") : u(x) = c G C"^} consisting of 
constant vector-valued functions. From (2.2) it follows that 

l|Dz||i^(0) ^ Cr^C2||z|||^(o), zeZ. (9.2) 

Due to compactness of the embedding H^{0;C^) C L2(C;C"), inequality 
of the form (9.2) can be satisfied only on a finitedimensional subspace in 
H^{0;C"'). Consequently, Z is finitedimensional. 

Obviously, the finitedimensional space (9.1) is also a (closed) subspace of 
L2(0;C"). We denote n{0) := L2(C>;C")e^, i?i(C>;C") = H\0;C'')n 
7i{0). In other words, 

ifi(0:C") = {u G Fi(0;C") : (u,z)i,(o) = 0, Vz G Z}. (9.3) 

Clearly, H]_{0;C"-) is a (closed) subspace of H^{0;C''). 

Proposition 9.1. The form ||6(D)u||^2(0) defines a norm in the space (9.3) 

equivalent to the standard H^-norm. 

Proof. The estimate ||6(D)u|||^(.q^ ^ aid||u||^ij^0^ obviously follows from 
(1.2) and (1.5) and is valid for all functions u G H^{0; C"). 

Let us prove the opposite estimate: there exists a constant Ci > such 
that 

l|u|||i(0)^Ci||6(D)u||i^(0), uGi/i(0;C"). (9.4) 

Assume the opposite. Suppose that for any A; G N there exists a function 
Ufe G H]_{0;C") such that 

l|ufcll|i{0) > ^IIKD)ufc||2^(o). 

Then for the sequence = ||ufc||^l^^jUfc we have ||vfc||j:/i(0) = 1 and 

1 1 6(D) Vfe 11/^2(0) < — ^ 0, as A; ^ oo. Due to compactness of the embed- 
ding ff^(0;C") C L2(0;C") there exists a subsequence {v^^.} converging 
in L2(C';C"). Applying the inequality (2.2) with u = v^. — v^, and using 
convergence of v^^ in L2{0) and convergence of 6(D) v/;^ in L2{0), we con- 
clude that { Vfc^ } is a fundamental sequence in H^{0;C'^). Consequently, it 
tends to some function v* G H]_{0;C'^). Then ||v*||jji(o) = 1, and, from 
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the other side, 5(D)v* = 0, i. e., v* G Z. Since Z n ifi(0;C") = {0}, we 
arrive at a contradiction. • 

Let An^e be the operator corresponding to the form (2.3). Ob- 
viously, Kcr^AT^e = Z. Therefore, the orthogonal decomposition 
L2(P\ C") = Z ® 7i{0) reduces the operator AN,e- Denote by l3N,e the self- 
adjoint operator in 7{(0) which is the part of AN,e in T-L{0). Equivalently, 
the operator Bn,^ is the selfadjoint operator in T-L{0) generated by the qua- 
dratic form 

6Ar,.[u, u] = / (/(x)6(D)u, 6(D)u)dx, u G Hi{0- C"). 
Jo 

This form is closed and, by (9.4), positive definite. We have 

ll5~iZl^rlu|l|i(0) < feiv,e[u,u] ^ ||5||Loedai||u||^i(0), u G Hi{0-C). 

(9.5) 

Our goal in this section is to approximate the inverse operator for 
small £. In terms of solutions, we study the behavior of the generalized 
solution Ue G H]_{0;C^) of the Neumann problem 

6(D)*/(x)6(D)u,(x) = F(x), x G O; 

d^Vi,\Qo = 0; (u„ z)i,(o) = 0, Vz G Z, ^ ' ^ 

where F G H(0). Then = B^^^F. 

By definition, we say that Ug G H]^{0; C") is the generalized solution of 
the problem (9.6) if it satisfies the identity 

/ (/6(D)u„6(D)r?)dx= / {F,rj)dy:, Vr? G /fi(0; C"). (9.7) 
Jo Jo 

By (9.5), the form tTv.efu, v] can be viewed as an inner product in the space 
H]_{0;C'^). The right-hand side of (9.7) is an antilinear continuous func- 
tional of r7 G H]_{0; C"). By the Riss theorem, the solution exists and is 
unique. It satisfies the estimate 

\Mh\o) ^ \\9~^\\l^Ci\\F\\l^^o)- (9.8) 

9.2. The homogenized problem. Let A^^^ be the operator gener- 
ated by the form (2.27). Obviously, Kcr^^ = Z. The decomposition 
L2iO; C") = Z® n{0) reduces the operator A%. Denote by B% the selfad- 
joint operator in H(0) which is the part of A% in T-L{0). In other words, 
is the selfadjoint operator in T-L{0) generated by the quadratic form 

6^[u,u] = /" (5°6(D)u,6(D)u)dx, uGi?i(0;C"). (9.9) 

Jo 

By (1.15), the form (9.9) satisfies twosided estimates of the form (9.5) with 
the middle part replaced by the form (9.9); the constants remain the same. 
Then the operator is positive definite: 

^N^Wg-^WllC^^Iuioy (9.10) 
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Let uq E H]_[0; C") be the generalized solution of the Neumann problem 

6(D)V6(D)uo(x)=F(x), x G O; 

d'luoldo = 0; (uo,z)i2(C') =0, Vz G Z, 

where F G n{0). Then uq = (.S^)"^F. 

The solution uq is defined similarly to (9.7). It exists, is unique, and 
satisfies the estimate 

l|uo|ki(o) ^ \\9~^\\lMMl2{0)- (9-12) 
By the theorem about regularity of solutions, we have Uq G H^{0;C'^) and 

WMhHO) ^ C\\F\\l^^o)- (9.13) 

Note that the estimate (9.13) can be derived from (2.33) and the fact that 
Uq is the solution of the problem (2.28) with the right-hand side F + Auq. 
Then, by (2.33) and (9.12), 

I|uo||h2(o) ^ C'aIIF + Auo||i2(c)) ^ (Ca + A||5~^||j:^Ci)||F||i2(o)- 
Fixing A = Ao = C2||s'~^||^^ + 1, we obtain (9.13) with the constant C = 
C*a„ + Ao||<7"^I|looCi. 

9.3. Approximation of solutions in L2(0;C"). Next theorem is an 
analog of Theorem 4.1 for the case where A = 0. 

Theorem 9.2. Suppose that O C is a bounded domain with the bound- 
ary of class C^'^. Suppose that the matrix g{'x.) and DO b(D) satisfy the 
assumptions of Subsection 1.2. Suppose that Condition 2.1 is satisfied. Let 
Us he the solution of the problem (9.6) and let Uq be the solution of the prob- 
lem (9.11) with F G HiO). Then there exists a number £i G (0, 1] depending 
on the domain O and the lattice T such that 

l|Ue - U0||l2(O;C") ^ Coe||F||i2(0;C")> < £ ^ £i, (9.14) 

or, in operator terms, 

The constant Cq depends on m, d, oq, ai, \\g\\LocJ lb ""^IliooJ the parameters of 
the lattice T, the constants Ci and C2 from the inequality (2.2), the constant 
Ci from (9.4), and the domain O. 

Proof. We apply Theorem 4.1. Let A be subject to the restriction (2.6). 
Since Bn^s is the part of AN,e in 'H{0) and is the part of in the same 
subspace, it follows from (4.2) that 

\\{BN,e + \ir^ -{Blj + \ir^\\Hio)^H{o)^CQWe, 0<e^£i. (9.15) 

Obviously, the solution of the problem (9.6) also satisfies the equa- 
tion {BN,e + A/)u£ = Fg, where F^ = F -|- Au^ G 7^(0). Consequently, 
Us = {BN,e + a/)~^F£. By (9.8), we have 

\\Fs\\l,{0) ^ \\nL,{0) + MMl,{0) ^ (1 + A||5~'||l^Ci)||F|U,(o). (9.16) 
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We put uo,e := {B% + A/)-1F,. By (9.15), 

l|ue - uo,£||l2(0) ^ Co(A)e||Fe||i2(0)> 0<e^ei. (9.17) 

From the other side, the solution uq of the problem (9.11) also satisfies 
the equation {B% + A7)uo = Fq, where Fq = F + Auq G H(0). Therefore, 
uo = {B% + A/)-iFo. We have: 

uo,e - Uo = {Bl + - Fo) = A(H^ + A7)-^(ue - Uo). (9.18) 

Consequently, 

Ue - Uo = U£ - Uo,£ + \{B%[ + A/)~-^(u£ - Uq), 

whence 

(7 - \{B% + A7)-i) (u, - uo) = u, - uo,e. (9-19) 

By (9.10), 

\\\{B% + \I)-^\\Hio)-,H{0) ^ A (a+ ||5-'|IZlCr^)"' < 1, 

and so the operator 7 — A(H^ + A7)~^ is invertible, and the norm of the 
inverse operator satisfies the estimate 

II (7 - X{Bl + Xiyy' |k(o)^^(o) ^ 1 + A||5-^||l^Ci. (9.20) 

As a result, (9.19) imphes that 

u, - uo = (7 - X{B% + A7)-i)"' (u, - UQ,,). (9.21) 

Combining (9.16), (9.17), (9.20), and (9.21), we obtain 

\\ne-uo\\L,(o) ^ (l + A||5-i||L^Ci)2Co(A)e||F||i,(o), < e ^ ei. (9.22) 

Estimate (9.22) is proved for any A satisfying (2.6). We can fix A, for in- 
stance, putting A = Ao := 1 + C2||5~^||Z^ ■ Then we obtain (9.14) with the 
constant Co = (1 + Ao||i^~^ ||l^Ci)2Co(Ao). • 

9.4. Approximation of solutions in 7f^(0;C"). The corrector for the 
problem (9.6) is defined similarly to (4.5): 

KN{e) = RomSeb(D)Po{B%r\ (9.23) 

The operator (9.23) is a continuous mapping of 7^(0) into H^{0; C"). The 
first order approximation to the solution of the problem (9.6) is given by 

V, = iB%)-'F + eKNie)F = uq + eA^5,6(D)5o, (9.24) 

where uq = PqUq. Note that belongs to H^{0; C*) but (in general) does 
not belong to 77|(C>;C"). 

The following result is an analog of Theorem 4.2 in the case where A = 0. 

Theorem 9.3. Suppose that the assumptions of Theorem 9.2 are satisfied. 
Let be defined by (9.24). Then there exists a number ei G (0, 1] depending 
on the domain O and the lattice T such that 

\\Ue - Ve||//i(0;Cn) ^ C£l/^||F||^2(o.cn), < £ ^ £i, (9.25) 
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or, in operator terms, 

For the flux := g^b(D)ue we have 

Hp, - g'S,b(P)uo\\L,iO;Crr^) ^ CV/'||F||i^(0;C"), < £ ^ £1, (9.26) 

where ^(x) is the matrix (1.9). The constants C and C depend only on m, 

d, ao, «!, IbllLoo) lb~"^llioo) the parameters of the lattice T, the constants Ci 
and C2 from (2.2), the constant Ci from (9.4), and the domain O. 

Proof. We apply Theorem 4.2. Suppose that A satisfies condition (2.6). De- 
note by V the orthogonal projection of L2{0; C") onto the subspace 7i{0). 
Multiplying operators under the norm sign in (4.10) by the projection V 
from the right, we arrive at 

\\{BN,e + A/)-^ - {B% + A7)-i - eA'S,b{-D)PoiB% + Xir'\\niO)^HHO) 

<C(A)£^/2, 0<£^£i. 

(9.27) 

As in the proof of Theorem 9.2, we use that = {BN,e + AI)"^F£ and 
consider the function uq,, = {B% + XI^^Fe- By (9.27) and (9.16), 

||u, - uo,s - eA'Seb{'D)uo,e\\HHo) ^ C{X)e^/^\\F,\\L^^o) 

^ (9.28) 

^ C(A)(1 + Mlg-'UM^'^^WnMO), < £ ^ £i, 

where uo,£ = Pouo,e- 

Note that from (2.32) and (2.34), by multiplying operators under the 
norm sign by the projection V from the right, it follows that 

\\iB% + XI)-'\\nio)^mio)^c'^\ (9-29) 

\\{B% + XI)-'\\n(^a)^HHO)^Cx. (9.30) 
Combining (9.18), (9.14), and (9.29), we obtain 

||U0,£ - UolUi(O) = A||(H^ + XI)-^{Ue - Uo)lliji(O) 
^ Ac^^llUe - Uo||l2(0) ^ Ac^^Co£ II F 11^2(0), < £ ^ £l. 

Now we consider the function £A^5e6(D)(uo,e — Uo)- It is required to 

estimate its norm in H^{0;C'^). Since this function is defined on the whole 
W^, we will estimate its norm in H^{M.'^;C^). By Proposition 1.5 and (1.4), 

||£A^S'£6(D)(uo,£ - uo)||L2(Md) ^ £|^r^^^ll^llL2(n)"}^^l|uo,£ - uo||//i(]Rd). 

(9.32) 

Consider the derivatives 

dj (£A^S,6(D)(So,e - 5o)) = (5, A)^S,6(D)(So,e - So) 

+ e A' S,b{r>)dj {uo,e-uo). 



(9.31) 
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By Proposition 1.5 and (1.4), we have: 
d 

J2\\dj{sA'S,b{B){uo,e-UomUm<i) 

^ 2|f^|"^||DA||^2(f2)ai||uo,e - uoll^i 
+ 2£2|J^|-i||A||2^(f^)ai||uo,e-5o||^2 
Together with (9.32) this imphes that 

||£A^5,6(D)(uo,.-uo)|||i(o) ^ ^ 

(9 33) 

^ |J^|-iai(3£2||A|||^^^^+2||DA||i^(^))||uo,e-uo||l,.(Krf). 
From (9.14), (9.18), and (9.30) it follows that 

||U0,£ — Uo||/f2(]Rd-) ^ CollUoe — Uo||ij2/o) 

^ ' ' ^ ~ (9.34) 

Ue — Uo||i2(o) ^ £Cc)ACaCo||F||^2(C')) < £ ^ £i. 

Combining (9.33), (9.34), (1.10), and (1.11), we obtain 

||£A^Se6(D)(uo,e - Uo)IUi(o) ^ C'A£||F||i^(o), < £ ^ £i, (9.35) 

whereC, = mV2aV2a-V2||^||V^||^-i||V^(2 + 3(2ro)-2)^/^CoAC,Co. 
As a result, relations (9.28), (9.31), and (9.35) imply that 

llus-VellHi(o) ^e'/'(c(A)(l + A|b-i||L^Ci) + Ac-iCo + C'A) ||F|U,(o) 

(9.36) 

for < £ ^ £i. Estimate (9.36) is proved for any A satisfying (2.6). Putting 
A = Ao := 1 + ^2115"^ 11^^ , we obtain estimate (9.25) with the constant 

C = C(Ao)(l + Ao|b-iL..?i) + AoCao'Co + Cx,. 

It remains to prove (9.26). The arguments are similar to the proof of 
(4.11) (see Subsection 5.4). By (9.25), (1.2), and (1.5), we have 

||p,-/6(D)v,||i,(0) ^ ||5|U^(dai)V2c£V2||F||^^((p), < £ ^ £i. (9.37) 

Similarly to (5.40), 

||/6(D)v, - rSeb(D)^\\L,iO) ^ Cee\\Bb{B)uo\\L,iR^y (9.38) 

By (1.4) and (9.13), 

||D6(D)uo||i2(M'') ^ "l''^l|uo||/f2( 



1/2 1/2 - (^•39) 

^ a/ CoWuqWh^o) ^ Oil CoC\\F\\l^(^o)- 

As a result, (9.37)-(9.39) imply (9.26) with the constant C' = 
\\g\\LAdai)^/^C + Ceal^^CoC. • 

9.5. Results in the case where A G Lqo. Under Condition 1.10, instead 
of the corrector (9.23), it is possible to use a simpler corrector 

<(£) = [A^]bmB%r\ 
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which is a continuous mapping of 'H{0) into i^^(0;C'^). Instead of (9.24), 
we consider another approximation of the solution of the problem (9.6): 

V, = {B%)-^¥ + eK%{e)¥ = uq + £A"6(D)uo. (9.40) 

Next result is an analog of Theorem 7.1 in the case where A = 0. 

Theorem 9.4. Suppose that the assumptions of Theorem 9.2 and Condition 
1.10 are satisfied. Let Vg he defined by (9.40). Then there exists a number 
£i G (0, 1] depending on the domain O and the lattice T such that 

l|U£ - ^e\\m{0;C^) ^ CV/2||F||^2(0;C"), < £ ^ £i, (9.41) 

or, in operator terms, 

For the flux Pe = 5^6(D)ue we have 

lIPe - 5'KD)uo||l2(0;C-) < CV/2||F||i,(o.cn), < £ ^ £i, (9.42) 

where ^(x) is the matrix (1.9). The constants C and C depend on m, d, oq, 
Q^i! Ilfflkoo) lb~"^IUoo) parameters of the lattice T, the constants Ci and 
C2 from the inequality (2.2), the constant Ci from (9.4), the norm ||A||i;,^, 
and the domain O. 

Proof. The proof is similar to that of Theorem 7.1. It is required to estimate 
the iJ^-norm of the function Ve — = eA^{I — Se)b{'D)uo. 
From (1.4) and (9.13) it follows that 

ll^(D)uo||i,2(Kd) ^ a}''^||uo||ii-i(Kd) ^ ay^Co\\uo\\H^(^o) < «i^^C'o^I|F||l2(C>)- 

(9.43) 

Estimates (7.6) and (7.7) can be applied. Together with (9.39) and (9.43) 
they yield 

l|Ve - Veki(o) = lkA^(/ - Se)b(B)uo\\HiiO) ^ C^124nL,{0), (9-44) 

where C?2 = a'/^CoC (2(3 + 2P2)^/^\\A\\l^ + (2/3i)i/Vi) . 

Now (9.25) and (9.44) imply (9.41) with the constant C = C + C^^- 
It remains to check (9.42). From (9.41), (1.2), and (1.5) it follows that 

||p,-/6(D)v,|U,(o) ^ £V2di/2„V2||<^||^^c||F|U,(o), < e £1. (9.45) 

Similarly to (7.10) and (7.11), taking (9.13) into account, we have 

||/6(D)ve -5^^6(D)uo||l2(ci) ^ £dai\\g\\L^\\A\\Lj\uo\\H2(o) 

^ (9.46) 
^ edai\\g\\Lj\A\\L^C\\F\\L^^o)- 

Relations (9.45) and (9.46) imply (9.42) with the constant C' = 
d'/'c,y'\\g\\Lj + dai\\g\\U\M\L^C. . 

9.6. Special cases. Now we distinguish special cases. Next statement 
follows from Theorem 9.3 and Proposition 1.2. 
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Proposition 9.5. Let Ug be the solution of the problem (9.6), and let uq 
be the solution of the problem (9.11). If =9, i- e., relations (1.13) are 
satisfied, then A = 0, Kn{£) = 0, and we have 

l|U£ — Uo||iji(C);C") ^ C£-''/^||F||i2(0;C")) < £ ^ £i. 

Applying Propositions 1.3, 1.13 and the statement of Theorem 9.4 con- 
cerning the flux, wc obtain the following statement. 

Proposition 9.6. Let be the solution of the problem (9.6), and let Uq 
be the solution of the problem (9.11). Let p^ = y'^6(D)u£. If = g, i. e., 
relations (1-14) are satisfied, then 

||Pe ~ /6(D)uo||i,(o.c-) CV/2||F||i^(o.cn), < e < £i. 

9.7. Approximation of solutions in a strictly interior subdomain. 

As in §8, it is possible to obtain a sharp order error estimate for approxima- 
tion of the solution in H^{0'; C"), where O' is a strictly interior subdomain 
of the domain O. 

Theorem 9.7. Suppose that the assumptions of Theorem 9.3 are sat- 
isfied. Let O' he a strictly interior subdomain of the domain O. Let 
5 := dist {O' ; dO} . Then there exists a number Si G (0,1] depending on 
the domain O and the lattice T such that 

- V£||ii-i(o';C") ^ ^5e||F||L2(0;C"), < £ ^ £i, (9.47) 
or, in operator terms, 

ll^ivje - {^nT^ - e-'^Jv(e)|k(0)->i?i(0';C") ^ ^i^, < £ ^ £i. 

For the flux p^ = 3^6(D)u£ we have 

\\Pe - WSeb{Ti)\io\\L^(^o';C^) ^ '^'6^\\P\\l2(0;C-'), < £ ^ £i, (9.48) 

where ^(x) is the matrix (1.9). The constants in estimates are given by 
£5 = ^\6~^ + ^2, ^'s = ^iS~^ + ^2, where £i, ^2, ^i, ^'2 depend on m, d, ao, 

Q^i) llffll-Loo) lb~"^llioo) parameters of the lattice T, the constants Ci and 
C2 from, the inequality (2.2), the constant Ci from (9.4), and the domain O. 

Proof. We apply Theorem 8.1. Assume that A satisfies (2.6). Multiplying 
operators under the norm sign in (8.2) by the projection V from the right, 
we arrive at 

\\{BN,e + A/)-^ - {B% + A/)-i - Ek'S,h{Ti)Po{Bl + \ir^\\H{0)^m{0') 

^ £5(A)£, < £ ^ £1. 

(9.49) 

Next, we proceed like in the proof of Theorem 9.3. We write as Ue = 
(BN,e + Ai")~^F£ and consider uo,e = + Xl^^^e- By (9.16) and (9.49), 

||u£ - uo,£ - £A^S'£6(D)uo,£||ii-i(ci') ^ (2:5(A)£||F£ 11^2(0) 

~ (9.50) 

^C5(A)(l + A||5-l||L^Ci)£||F||i2(0), 0<£^£i. 
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Now relations (9.31), (9.35), and (9.50) imply that 

Ilue - Mm{0') ^ £ (^6(A)(1 + M\9~^\\l^Ci) + Xc^^Co + Cx) \\F\\l2{0) 

(9.51) 

for < e ^ £i. The inequality (9.51) is valid for any A satisfying (2.6). 
Putting A = Ao = 1 + C2||5~''^||£^, we obtain (9.47) with the constant = 
£i(5~^ + £2, where 

£1 = £i(Ao)(l + Aoll^-^L^Ci), (9.52) 

£2 = £2(Ao)(l + XoWg-^hM + Aoc^o'Co + Cxo- 

It remains to prove (9.48). From (9.47), (1.2), and (1.5) it follows that 

lIPe - /6(D)v,|U,(o,) ^ \\g\\L^d^/^ay^ls£\\nL,{0), < £ ^ £1. (9.53) 
Combining (9.53), (9.38), and (9.39), we arrive at 

Hp, - rSeHB)uo\\L2(0') ^ {\\g\\L^d^/^a/^ls + CQa]/^CoC) e||F||i,(o) 

for < £ ^ £1. This yields (9.48) with the constant £'5 = £i5-^ + £'2, where 

= ll5l|Loof^'/'«}^'^i, l'2 = \\g\\L^d^'''a^'''^2 + C^a]'''CoC. . (9.54) 

Under Condition 1.10, we obtain the following result. 

Theorem 9.8. Suppose that the assumptions of Theorem 9.4 are satis- 
fied. Let O' he a strictly interior suhdomain of the domain O, and let 
5 := dist{0';90}. Then there exists a number £1 G (0,1] depending on 
the domain O and the lattice V such that 

l|U£ - ^e\\m{0';C") ^ ^5e||F||L2(0;Cn), < £ ^ £1, (9.55) 

or, in operator terms, 

\\^N]e ~ (^Af)""^ ~ £K%{e)\\'^(o)^m{0';C^) ^ ^<5^> < £ ^ £1. 

For the flux Pe = 5('^6(D)ue we have 

||Pe - F^(D)uo||l2(C)';C-) ^ ^5e||F||L2(0;C-), < £ ^ £l, (9.56) 

where g{x) is the matrix (1.9). The constants in estimates are given by 
£5 = £i(5~^ + £2, £5 = £'i<5^^ + ^2' where £1, £2, £'1, £2 depend on m, d, ao, 
o^i) llfflUoo) ll5~^IUoo) parameters of the lattice F, the constants Ci and 
C2 from the inequality (2.2), the constant Ci from (9.4), the norm ||A||/,^, 

and the domain O. 

Proof. We apply the proof of Theorem 9.4. Under Condition 1.10, (9.44) is 
true. Together with (9.47) it implies (9.55) with the constant £5 = £i(5^^ + 
£2, where £1 is defined by (9.52) and £2 = £2 + C^2- 

Let us check (9.56). Prom (9.55), (1.2), and (1.5) it follows that 

Hp, - /6(D)v,|U,(o,) ^ \\g\\L^d^l''a\'^'^5e\\F\\L^^o), < £ ^ £1. (9.57) 
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Under Condition 1.10, (9.46) is true. Now relations (9.46) and (9.57) implies 
(9.56) with the constant C'^ = ^[6-^ + €'2, where is defined by (9.54) and 
^2 = \\9\\L^d'^^a/^C2 + dai\\g\\Lj\A\\L^d. . 

What was said in Remark 8.3 concerns also Theorems 9.7 and 9.8. 
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